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1.0  SCOPE 


1.1  Identification 

This  Software  Design  Document  (SDD),  prepared  for  the  Oceanographic  and  Atmospheric 
Master  Library  (OAML),  provides  detailed  information  on  the  nearshore  wave  and  current 
forecasting  software  titled  SURF  3.1.  This  model  equips  users  with  an  automated  method  for 
determining  surf  conditions  and  related  environmental  parameters.  SURF  3.1  produces  a  standard 
surf  forecast  and  a  Modified  Surf  Index  (MSI)  number,  which  are  Navy  requirements  for  littoral 
operations  and  amphibious  landings  (see  Joint  Surf  Manual).  The  first  operational  Navy  surf 
forecasting  computer  model  was  developed  for  the  Fleet  in  1988  (see  Earle,  1988)  to  supplement  the 
manual  and  visual  techniques  developed  in  the  1950's.  The  manual  procedures  are  subjective  and 
do  not  adequately  consider  shallow  water  effects  such  as  wave  shoaling  and  refraction.  This  version 
of  SURF  3.1  is  a  modem  redesigned  application,  which  uses  state  of  the  art  technology  in  operational 
real-time  surf  zone  forecasting. 

1.2  Document  Overview 

This  OAML  SDD  describes  the  design,  structure,  and  scientific  aspects  of  the  Computer 
Software  Configuration  Item  (CSCI)  titled  SURF  3.1.  This  document  provides  a  detailed  summary 
of  all  Computer  Software  Units  (CSU)  or  subroutines,  input  file  formats,  output  file  formats,  and 
user-specified  options.  The  SDD  is  divided  into  three  sections;  the  Preliminary  Design,  the 
Architectural  Design,  and  the  Detailed  Design. 

The  Preliminary  Design  section  describes  the  scientific  aspects  of  SURF  3.1  including  a  brief 
description  of  the  mathematical  formulation  and  theory  behind  the  model.  The  Architectural  Design 
section  outlines  the  structural  design  of  SURF  3.1  with  a  graphical  representation  of  the  CSU  calling 
sequence.  The  Detailed  Design  section  identifies  and  summarizes  the  operation  of  each  CSU 
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including  detailed  listings  of  input  variables,  output  variables,  local  variables,  calling  routines,  and 
called  routines  and/or  called  functions. 
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3.0  PRELIMINARY  DESIGN  OF  SURF  3.1 


3.1  SURF  3.1  CSCI  Overview 

SURF  3.1  is  a  parametric  one-dimensional  model  based  largely  on  the  work  of  Thornton  and 
Guza  (1983,  1986).  Thornton  and  Guza  developed  several  models  for  random  wave  processes 
including  a  wave  height  transformation  model  and  a  longshore  current  model.  These  models  contain 
both  numerical  and  analytical  solutions,  which  provide  cross-shore  distributions  of  various 
parameters  such  as  wave  height,  longshore  current  velocity,  and  wave  length.  However,  because 
SURF  3.1  is  one-dimensional,  certain  approximations  are  made:  (1)  straight  and  parallel  bottom 
contours,  (2)  depth-uniform  currents,  (3)  wave  heights  are  Rayleigh  distributed,  (4)  linear  wave 
theory  is  applicable,  and  (5)  directional  wave  spectra  are  narrow-banded  in  frequency  and  direction. 

The  model  is  designed  to  operate  in  a  variety  of  modes  to  provide  both  military  and  civilian 
users  with  local  surf  and  current  forecasts.  SURF  3.1  requires  three  distinct  pieces  of  information 
to  perform  calculations:  (1)  a  depth  profile,  (2)  a  directional  wave  spectrum,  and  (3)  wave  refraction 
information.  Each  of  these  required  data  sources  can  be  accessed  externally  or  generated  internally. 
This  design  allows  for  maximum  flexibility  when  using  SURF  3. 1  to  generate  local  forecasts  where 
input  data  may  or  may  not  be  available.  The  details  of  these  input  formats  are  described  in  section 
6.0.  The  following  subsections  outline  the  scientific  principles  behind  SURF  3.1  and  the  inherent 
fundamental  hydrodynamic  calculations  contained  in  the  model. 

3.1.1  Wave  and  Roller  Energy  Models 

As  waves  approach  the  coast,  the  frictional  effect  of  the  sea  floor  on  the  organized  orbital 
motion  of  water  particles  within  a  wave  causes  the  wave  to  break  or  spill.  The  flows  of  spilling 
breakers  can  be  separated  into  two  layers,  an  upper  layer  of  turbulent  energy,  which  rides  over  a 
lower  layer  of  energy  that  maintains  an  organized  oscillatory  wave  motion.  The  region  of  turbulent 
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water  above  the  wave  is  termed  a  surface  roller.  The  original  idea  of  such  a  two-layer  system  was 
introduced  by  Longuet-Higgins  and  Turner  (1974)  (see  also  Svendsen  (1984a,b)).  SURF  3.1 
incorporates  the  model  of  Lippman  et  al  (1995,  1996),  which  produces  results  consistent  with 
measurements  from  both  a  planar  and  a  barred  beach.  The  energy  associated  with  each  region  of 
interest  is  utilized  to  shoal  the  incoming  waves  and  drive  the  longshore  current.  The  energy  per  unit 
surface  area  in  a  wave  is  calculated  as: 


Ew  =  ~pgHL 

O 

where  p  is  water  density  and  g  is  the  acceleration  due  to  gravity.  Hnm  is  the  root-mean-square  wave 
height.  The  energy  per  unit  area  associated  with  a  roller  is  given  as: 


Er  = 


Hi 

h  tan  u 


where  c  is  the  phase  speed  of  the  wave, /is  the  zero  crossing  frequency,  Hb  is  the  height  of  the  wave 
at  breaking,  h  is  water  depth,  and  a  is  the  angle  the  roller  makes  with  the  body  of  the  wave.  A 
default  value  of  5  degrees  is  used  for  the  roller  angle  in  SURF  3. 1 . 

3.1.2  Energy  Dissipation  in  the  Surf  Zone 

As  a  wave  propagates  across  the  surf  zone,  its  energy  is  dissipated  due  to  bottom  friction, 
wave  breaking,  turbulence,  and  wave-current  interaction.  A  generic  formulation  of  this  energy 
dissipation  is  given  by  the  energy  flux  equation: 

d(E„cg  cos  0)  _ 
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where  Ew  is  the  wave  energy,  cs  is  the  wave  group  velocity  and  6  is  the  wave  direction  relative  to 
shore  normal  (x  positive  offshore).  The  Right  Hand  Side  (RHS)  of  the  above,  equation,  <£p%  is  the 
ensemble  averaged  dissipation  function.  Thornton  and  Guza  (1983)  modeled  this  dissipation 
function  as: 


4  h 

where  B  is  an  empirical  coefficient,  and  pb(H)  is  the  probability  distribution  for  breaking  waves 
described  by: 


pb(H)  =  W(H)p(H) 

where  p(H)  is  a  Rayleigh  Distribution  of  wave  heights  and  W(H)  is  a  weighting  function  resulting 
in  a  weighted  Rayleigh  distribution.  Several  weighting  functions  W(H)  have  been  constructed  by 
various  authors,  the  weighting  function  applied  in  SURF  3.1  developed  by  Thornton  and  Guza 
(1986)  is  given  as: 


l-e 

v  J 

where  y  is  an  empirical  factor  determined  from  field  data  to  be  0.42,  h  is  the  water  depth  and  H  is 
the  wave  height.  If  wave  roller  energy  is  considered  in  the  model,  the  modified  energy  flux  equation 
is  given  as: 


w{h)= 


H. 


-i4 


y  h 


JL 

yh 


d(  E w  cg  cos  6)  d(ErCCOs6)_  = 


dx 


dx 


=  -<Sr> 
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and  the  dissipation  becomes  a  function  of  the  roller  term. 


<£>  =  ^pgf^coscj  \ H3 Pb(H)dH 
4  h 

The  above  equation  is  solved  using  a  numerical  forward  stepping  and  convergence  scheme  to 
determine  wave  and  roller  energy  along  with  Hrms  values  at  each  point. 


3.1.3  Longshore  Current  Calculations 

When  waves  enter  the  surf  zone  at  an  angle,  the  shore-parallel  component  of  momentum 
inherent  to  wave  motion  drives  a  current  along  the  shore.  This  longshore  current  can  be  a  significant 
force  inside  the  surf  zone.  Calculation  of  the  current  velocity  is  based  on  radiation  stress  theory  (see 
Longuet-Higgins,  1970a,  1970b).  A  general  form  of  the  longshore  momentum  equation  is: 


( 


Ty+P  — 

dx 


/  dV 


-<r?>  +  r  ;  =  0 


where  p  is  the  water  density,  h  is  the  water  depth,  and  V  is  the  longshore  current.  The  first  term  on 
the  left  hand  side  is  the  radiation  stress  in  the  along  shore  direction  exerted  by  waves  on  the  water 
given  by: 


sin# 


where  eb  is  the  dissipation  function  defined  in  the  previous  section,  c  is  wave  phase  speed,  and  G  is 
the  angle  of  wave  approach  with  respect  to  x.  The  second  term  is  the  horizontal  mixing.  The 
horizontal  eddy  viscosity  p  is  modeled  after  Battjes  (1975). 
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in  which  M  is  an  empirical  constant  equal  to  2.  The  third  term  is  the  mean  stress  due  to  bottom 
friction  given  by: 


Ty~  P  CfUV 

where  c/is  the  bottom  friction  coefficient,  u  is  the  magnitude  of  the  near-bottom  horizontal  wave 
orbital  velocity,  and  V  is  the  longshore  current.  Linear  wave  theory  defines  the  near-bottom  wave- 
induced  orbital  velocity  as: 


xH_ 

U  T  sinh  (kh ) 

where  H  is  the  wave  height,  T  is  the  wave  period  and  k  is  the  wave  number  which  can  be  calculated 
using  the  dispersion  relation: 

cr2  =  g  k  tanh  (k  h) 

where  cr  is  the  radian  wave  frequency  and  g  is  gravity.  The  longshore  current  equation  is  solved 
using  a  finite  difference  approach  after  wave  heights,  water  depths,  and  wave  dissipation  values  are 
calculated  at  each  cross-shore  grid  point  in  the  surf  zone. 

A  major  improvement  to  the  longshore  current  calculation  is  included  in  Surf  3. 1 .  Hsu  et  al. 
(2000)  showed  that  using  a  variable  bottom  friction  coefficient  in  the  longshore  current  model 
provides  more  realistic  distributions  of  longshore  current  velocities.  The  depth  dependent  bottom 
friction  coefficient  function  is  defined  as 


cf(x)  = 


0.003 


;  x> 


0.003 


w> 

V  J 


,x< 


X „ 


where  x  is  the  offshore  distance,  h  is  the  local  water  depth,  and  Xb  is  the  distance  from  the  shoreline 
to  the  location  where  ten  percent  of  the  waves  are  breaking.  It  should  be  noted  that  the  variable 
bottom  friction  function  reflects  the  shoreward  increase  in  friction  due  to  sediment  sorting  and 
compensates  for  the  lack  of  vertical  diffusivity  in  one-dimensional  models. 

3.1.4  Directional  Energy  Spectra 

SURF  3.1  allows  users  to  generate  surf  forecasts  using  two  different  directional  wave  energy 
spectrum  types.  The  user  can  choose  from  an  internally  generated  wave  spectrum  or  an  external 
wave  spectrum.  If  the  internally  generated  spectrum  is  selected,  a  modified  Pierson-Moskowitz 
(1964)  spectrum  is  calculated  based  on  sea  and  swell  conditions  defined  in  the  surf  model  input  file. 
A  detailed  description  of  the  external  wave  spectrum  format  is  available  in  section  6.0.  Users  can 
also  examine  the  shoaled  and  refracted  directional  wave  spectrum  at  specific  depths  by  using  options 
described  in  section  6.0. 

3.1.5  Differences  Between  Surf  3.0  and  Surf  3.1 

The  transition  from  Surf  3.0  to  Surf  3.1  includes  several  scientific  and  code  improvements. 
Surf  3.1  includes  a  new  longshore  current  model,  based  on  the  work  of  Hsu  et  al.  (2000),  which 
provides  improved  longshore  current  velocity  distributions.  Several  error  checking  routines  to 
examine  the  stability  and  usability  of  input  depth  profiles  are  used  in  Surf  3. 1 .  The  length  of  input 
file  names  was  increased  to  40  characters.  The  ability  to  output  a  shoaled  and  refracted  shallow 
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water  directional  wave  spectrum  is  now  available  to  expert  users.  A  description  of  the  output 
directional  wave  spectrum  file  is  available  in  section  6.0. 

The  input  file  for  the  new  model  is  simplified  and  streamlined.  The  number  of  lines  in  the 
Surf  3.1  input  file  has  been  reduced  and  the  method  for  modifying  user  options  has  changed.  All  of 
the  options  are  still  available,  but  normal  operation  of  Surf  3.1  no  longer  requires  the  user  to  specify 
every  option  flag.  To  simplify  the  input  file,  a  set  of  default  options  is  implicitly  included  in  the 
input  file  format.  The  changes  to  the  input  file  are  described  in  section  6.0. 
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4.0  SURF  3.1  CSCI  ARCHITECTURAL  DESIGN 


The  Architectural  Design  section  shows  the  overall  design  and  the  calling  sequence  for  all 
CSU’s  of  the  SURF  3.1  Model.  Each  CSU  is  shown  in  the  calling  sequence  with  the  associated  CSU 
related  to  each  specific  unit.  Figure  (1)  presents  the  path  in  which  each  CSU  is  called  and  all 
associated  CSU’s,  which  in  turn  are  called  from  the  parent  unit.  Specific  details  concerning  the 
criteria  for  each  CSU  being  called  are  defined  in  the  Section  5.0:  SURF  3.1  CSU  Detailed  Design. 
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Figure  1.  SURF  CSCI  Architectural  Design 


WEIGHTFN 


5.0  SURF  3.1  CSCI  DETAILED  DESIGN 


5.1  Program  SURF 
Program  Call: 

SURF ( ) 


Summary: 


The  SURF  routine  is  the  starting  point  for  executing  SURF  3.1.  The  routine  identifies  the 


input  type  and  controls  the  reading  of  data  and  user  selected  computation  options.  The  routine  calls 
the  main  wave  parameter  calculation  routines  and  controls  the  output  of  the  resulting  data. 

Input  Variables:  None. 

Output  Variables:  None. 

Local  Variables: 


alfa 

Real 

bravo 

Real 

chrlie 

Real 

dangle 

Real 

depname 

Char*40 

dsea 

Real 

dstart 

Real 

dswell 

Real 

dxyl  (points) 

Real 

echo 

Real 

ehsig 

Real 

esowm  (dirNum,freqNum) 

Real 

filedat 

Char*40 

filein 

Char*40 

file_out 

Char*40 

filetmp 

Char*40 

foxtrt 

Real 

fracname 

Char*40 

freq  (freqNum) 

Real 

freql  (freqNum) 

Real 

freq2  (freqNum) 

Real 

Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Angle  Between  Directional  Bins 
Depth  Profile  File  Name 
Input  Direction  for  Sea  Contribution 
Input  Starting  Depth 

Input  Swell  Direction  for  Internally  Generated 
Spectrum 

Corresponding  Depths  with  No  Tide 
Breaker  Angle 

Significant  Wave  Height  from  Directional 
Spectrum 

Directional  Wave  Spectrum 
Output  File  Name  *.dat 
Input  Filename 
Output  File  Name  *.out 
Temporary  File 

Longshore  Current  Speed  and  Direction 
Wave  Refraction  File  Name 
Input  Wave  Spectrum  Center  Frequencies 
Beginning  Frequency  Bin  Values 
Ending  Frequency  Bin  Values 
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gamma2 

Real 

golfl 

Real 

golf2 

Real 

gt_frq 

Integer 

hsea 

Real 

hswell 

Real 

iday 

Integer 

idirec 

Integer 

ifreq 

Integer 

igamma 

Integer 

ihour 

Integer 

ihtll 

Real 

ihtl2 

Real 

imin 

Integer 

imonth 

Integer 

iyear 

Integer 

j  gamma 

Integer 

line 

Char*  80 

Unstress 

Logical 

lndname 

Char*40 

nnn 

Integer 

pet  (4) 

Real 

period  (freqNum) 

Real 

psea 

Real 

pswell 

Real 

roller 

Logical 

selfst 

Char*l 

slope 

Real 

spectra 

Logical 

spefile 

Char*40 

surfy 

Logical 

tide 

Real 

wdir 

Real 

Beach  Orientation,  Compass  Heading  Directly 

Toward  Beach 

Number  of  Surf  Lines 

Surf  Zone  Width 

Spectrum  Type 

Input  Significant  Wave  Height  for  Sea 
Contribution  to  Pierson  Moskowitz  Spectrum 
Input  Significant  Wave  Height  for  Internally 
Generated  Spectrum 
Input  Day 

Number  of  Direction  Bins  in  the  Input  Spectrum 
Number  of  Frequency  Bands  in  the  Input 
Spectrum 

Beach  Orientation  Rotated  90°  from  Original 
Heading  Toward  Beach 
Input  Hour 

Wind  Speed  Coded  Surf  Forecast  Value 

Wind  Direction 

Input  Minute 

Input  Month 

Input  Year 

Temporary  Value  Set  to  Beach  Orientation 

Temporary  Character  Variable 

Longshore  Current  Solution  (Tme  or  False) 

Input  Landing  Zone  Name 

Number  of  Points  in  the  Input  Depth  Array 

Percent  of  Different  Breaker  Types 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

Period  Array  (1 /Frequency) 

Input  Wave  Period  for  Sea  Contribution  to 

Pierson  Moskowitz  Spectrum 

Input  Swell  Period  for  Internally  Generated 

Spectrum 

Roller  Usage  (True  or  False) 

Self  Start  Flag  (Yes  or  No) 

Bottom  Slope 

Does  Input  Spectra  Exist?  (Tme  or  False) 
Selected  Wave  Spectrum  File  Name 
Significant  Wave  Heights  Greater  than 
0.5  ft?  (Tme  or  False) 

Input  Tide  Level 

Input  Wind  Direction  Compass  Heading  Wind 
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Blows  from 

wspd 

Real 

Input  Wind  Speed 

xcoeff  (dirNum,  freqNum) 

Real 

Wave  Height  Refraction  Coefficients 

xdelt 

Real 

Surf  Zone  Output  Interval 

xdeltgr 

Real 

Self  Adjusting  Cross-Shore  Grid  Step 

xfrom  (dirNum) 

Real 

Direction  Array,  Direction  Wave  Energy  Comes 
From 

xtheta  (dirNum, freqNum) 

Real 

Angle  Refraction  Coefficients 

xx  1  (points) 

Real 

Adjusted  Cross-Shore  Distances  from 

Depth  Profile 

ydepth 

Char*  1 

Input  Depth  Profile  Used?  (Yes  or  No) 

ydetail 

Char*  1 

Detailed  Output?  (Yes  or  No) 

yrefrac 

Char*  1 

Is  Refraction  Considered  in  Analysis? 

(Yes  or  No) 

ystr 

Char*  1 

Is  Straight  Coast  Refraction  Used? 

(Yes  or  No) 

Subroutines  Called  from  SURF  ( ): 

ABORT 

BREAD 

CALCSURF 

DEPDRVR 

GENSPEC 

MDSRF1 

PRTOUT3 

READRFRC 

READSPEC 

REFRAC 

SRFSETUP 

STRFRAC 

SUMMARY 

SURFCAST 
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Figure  2.  Program  SURF  Flowchart 
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5.2  Subroutine  ABORT 


Subroutine  Call: 

ABORT ( ) 

Summary: 

Subroutine  ABORT  acts  as  the  single  program  termination  routine.  The  subroutine  handles 
normal  program  execution  and  error  interrupt.  ABORT  is  called  to  stop  the  execution  of  the 
program.  If  an  error  interrupt  calls  ABORT  the  error  message  is  generated  locally  in  the  calling 
routine. 

Input  Variables:  None. 

Output  Variables:  None. 

Local  Variables:  None. 

Subroutines  Called  from  ABORT  ( ):  None. 

ABORT  (  )  Called  from  Subroutines: 

CJNJDEP 

EQUILPRF 

MAIN_WAV 

MDSRF2 

NEWBRK 

NONLIN2 

PRTOUT1 

PRT_OUT2 

READRFRC 

READSPEC 

SRFSETUP 

SURF 


Page  18 


Figure  3.  Subroutine  ABORT  Flowchart 


Page  19 


5.3  Subroutine  B_DETAIL 
Subroutine  Call: 

BDETAIL  ( iyear,  imonth,  iday,  ihour,  imin) 

Summary: 

Subroutine  B  DETAIL  formats  and  writes  the  detailed  surf  model  data  output  to  the  output 
text  file.  The  file  name  is  generated  as  “Lout”,  where  the  is  replaced  with  the  prefix  of  the  input 
file  name. 

Input  Variables: 


iday 

Integer 

Input  Day 

ihour 

Integer 

Input  Hour 

imin 

Integer 

Input  Minute 

imonth 

Integer 

Input  Month 

iyear 

Integer 

Input  Year 

Output  Variables: 

None. 

Local  Variables: 

None. 

Subroutines  Called  from  B  DETAIL  ( ):  None. 

B  DETAIL  ( )  Called  from  Subroutines: 

CALCSURF 

Figure  4.  Subroutine  B  DETAIL  Flowchart 


START 


I 


Writing  Header 
Information  to  the 
Output  File 


RETURN 
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5.4  Subroutine  B  HEAD 


Subroutine  Call: 

BHEAD  (gt_ffq,  roller,  Un  stress) 

Summary: 

Subroutine  B  HEAD  writes  header  information  and  user  selected  model  options  to  the 


output  file. 

Input  Variables: 


gtfrq  Integer 

lin_stress  Logical 

roller  Logical 

Output  Variables:  None. 

Local  Variables:  None. 


Subroutines  Called  from  B  HEAD  ( ): 

B  HEAD  ( )  Called  from  Subroutines: 
SURF 

Figure  5.  Subroutine  B  HEAD  Flowchart 


Spectrum  Type 

Longshore  Current  Solution  (True  or  False) 
Roller  Option  Flag  (True  or  False) 


None. 


START 

T 


Writing  Selected 
r  Input  Options  to  they 
Output  File 


_L 

RETURN 
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5.5  Subroutine  BALANCEQ 


Subroutine  Call: 

BALANCEQ  (roller,  theta,  Cg,  rhs,  hrmsl,  dp,  meanfreq,  xk,  hrms2,  convg) 

Summary: 

Subroutine  BALANCEQ  computes  the  new  wave  height  value  at  the  next  onshore  grid  cell 
by  performing  an  iterative  solution  to  the  energy  equations. 

Input  Variables: 


Cg 

Real 

Wave  Group  Velocity 

dp 

Real 

Offshore  Water  Depth 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 

mean_frq 

Real 

Wave  Frequency 

rhs 

Real 

Right  Hand  Side  of  Energy  Balance 

Equation 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

theta 

Real 

Wave  Angle 

xk 

Real 

Wave  Number 

Output  Variables: 

convg 

Logical 

Convergence  Flag  (True  or  False) 

hrms2 

Real 

Significant  Wave  Height  at  next 

Onshore  Grid 

Local  Variables: 

avgh 

Real 

Average  Wave  Height 

check 

Real 

Convergence  Check 

done 

Logical 

Flag  indicating  End  of  Loop 

f3 

Real 

Function  which  Calculates  Total  Energy 

kount 

Integer 

Loop  Iteration  Counter 

lhs 

Real 

Left  Hand  Side  of  the  Energy  Equation 

limit 

Logical 

Flag  for  Comparison  of  the  Left  &  Right  Side 
of  the  Energy  Equation  (True  or  False) 

lowerh 

Real 

Lower  Limit  of  Wave  Height 

maxkount 

Integer 

Maximum  Number  of  Loop  Iterations  =1000 

oldavgh 

Real 

Previous  Average  Wave  Height  Value 

pet 

Real 

Convergence  Step  Value 

tol 

Real 

Convergence  Tolerance 
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Upper  Limit  of  the  Wave  Height 


upperh  Real 

Subroutines  Called  from  BALANCEQ  ( ):None 
Functions  Called  from  BALANCEQ  ( ): 

F3 

BALANCEQ  ( )  Called  from  Subroutines: 

MAEN_WAV 

SLF_STRT 
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5.6  Subroutine  C  FINE 


Subroutine  Call: 

CJFINE  (ndepth,  xxin,  zzin,  xdelt  gr,  nnn,  xxl,  dxyl) 


Summary: 

Subroutine  C_FTNE  linearly  interpolates  the  input  water  depths  and  offshore  distances  to  an 
evenly  spaced  grid.  The  internally  defined  grid  self-adjusts  to  maximize  spatial  resolution  without 


exceeding  the  array  dimensions. 

Input  Variables: 


ndepth  Integer 

xdelt_gr  Real 

xxin  (points)  Real 

zzin  (points)  Real 

Output  Variables: 

dxyl  (points)  Real 

nnn  Integer 

xxl  (points)  Real 


Local  Variables: 


dxl 

dx2 

dxx 

dzz 

mm 

mml 

mmm 

nn 

xlast 

xtemp 


Real 

Real 

Real 

Real 

Integer 

Integer 

Integer 

Integer 

Real 

Real 


Subroutines  Called  from  C_FINE  (  ): 


Number  of  Points  in  Input  Depth  Profile 
Self-Adjusting  Cross-Shore  Grid  Step 
Cross-Shore  Distances 
Corresponding  Depths 


Corresponding  Depths  without  Tide 
Number  of  Points  in  the  Input  Depth  Array 
Adjusted  Cross-Shore  Distances  from 
Depth  Profile 


Temporary  Variable  Used  in  Calculation  of 

Next  Grid  Point  Distance 

Temporary  Variable  Used  in  Calculation  of 

Next  Grid  Point  Distance 

Distance  Quotient 

Difference  Between  Depth  and 

Previous  Depth 

Counter  Variable 

Counter  Variable 

Counter  Variable 

Counter  Variable 

Last  Distance  Offshore  from  Input  Profile 
Temporary  Variable  for  Cross-Shore  Values 

None. 
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C_FINE  ( )  Called  from  Subroutines: 

CREGRID 

Figure  7.  Subroutine  C  FINE  Flowchart 
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5.7  Subroutine  C  GAMMA 


Subroutine  Cali: 

CGAMMA  (gamrnaZ,  igamma) 

Summary: 

Subroutine  C  GAMMA  rotates  the  beach  orientation  data  read  from  the  input  file.  The  user 
defines  the  beach  orientation  as  the  compass  heading  of  a  boat  traveling  directly  toward  the  shore 
on  a  perpendicular  line  to  the  coast.  The  value  is  then  rotated  to  reflect  the  orientation  of  the  local 
coastline  with  respect  to  magnetic  north. 

Input  Variables: 

gamma2  Real  Beach  Orientation,  Heading  Directly 

Toward  Beach 

Output  Variables: 

igamma  Integer  Rotated  Beach  Orientation 

Local  Variables: 

gammatp  Real  Temporary  Variable  Used  in  Calculation 

mtemp  Integer  Temporary  Variable  in  Calculation 

Subroutines  Called  from  C_GAMMA  ( ):  None. 

C  GAMMA  ( )  Called  from  Subroutines: 

DEPDRVR 

Figure  8.  Subroutine  C  GAMMA  Flowchart 
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5.8  Subroutine  CINDEP 
Subroutine  Call: 

C_IN_DEP  (depname,  dstart,  xdeltgr,  nnn,  xxl,  dxyl) 

Summary: 

Subroutine  C  IN  DEP  reads  the  depth  profile  and  header  information  contained  in  the  input 
data  file.  The  routine  identifies  the  units  of  measurement  used  to  construct  the  depth  profile  and 
checks  the  order  of  the  offshore  distances.  If  the  data  is  misaligned,  the  subroutine  will  sort  and 
reorder  the  depths  and  offshore  distances. 


Input  Variables: 


depname 

Char*40 

Depth  Profile  File  Name 

dstart 

Real 

Input  Starting  Depth 

xdelt_gr 

Real 

Self  Adjusting  Cross-Shore  Grid  Step 

Output  Variables: 

dxyl  (points) 

Real 

Corresponding  Depths  without  Tide 

nnn 

Integer 

Number  of  Points  in  the  Input  Depth  Array 

xxl  (points) 

Real 

Adjusted  Cross-Shore  Distances  from  the  Depth 
Profile 

Local  Variables: 

al 

Real 

Temporary  Variable 

a2 

Real 

Temporary'  Variable 

adum 

Char*  80 

Temporary'  Variable,  Character  String  in  Input 
Field 

dcall 

Real 

Conversion  Factor  for  Distance  Offshore, 
Convert  to  Meters 

dcal2 

Real 

Conversion  Factor  for  Depths  Offshore,  Convert 
to  Meters 

dx 

Real 

Temporary  Variable  for  Distance  Offshore  from 
Input  File 

dz 

Real 

Temporary  Variable  for  Depths 

I 

Integer 

Loop  Variables 
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icall 

Integer 

ical2 

Integer 

instat 

Integer 

j 

Integer 

k 

Integer 

line 

Integer 

loop 

Integer 

ndepth 

Integer 

xxin  (points) 

Real 

zzin  (points) 

Real 

Subroutines  Called  from  C_IN_DEP  ( ): 

ABORT 

C_UN 

CREGRID 

C_IN_DEP  ( )  Called  from  Subroutines: 

DEPDRVR 


Input  from  Depth  File, 

Units  of  Distance  Offshore 

1  =  Feet 

2  =  Meters 

3  =  Yards 

Depth  Units  Input  from  Depth  File 

1  =  Feet 

2  =  Meters 

3  =  Fathoms 
File  Open  Status 
Loop  Variables 

Temporary  Variable  for  Number  of  Points 
Counter  for  the  Number  of  Lines  in  the  Input 
Depth  Profile 
Loop  Counter 

Number  of  Points  in  Input  Depth  Profile 
Cross-Shore  Distances 
Corresponding  Depths 
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Figure  9.  Subroutine  C  IN  DEP  Flowchart 


!. 
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5.9  Subroutine  C  REGRID 


Subroutine  Call: 

CREGRID  (ndepth,  xxin,  zzin,  xdelt  gr,  nnn,  xxl,  dxyl) 

Summary: 

Subroutine  C  REGRID  examines  the  cross-shore  step  size  (Ax)  of  the  input  depth  profile 
and  selects  a  new  step  size  to  optimize  the  depth  and  cross-shore  distance  arrays.  The  step  size  is 
automatically  adjusted  and  the  arrays  are  constructed  so  the  length  does  not  exceed  the  dimension 
of  the  array. 


Input  Variables: 


ndepth 

Integer 

Number  of  Points  in  Depth  Profile 

xdelt_gr 

Real 

Self  Adjusting  Cross-Shore  Grid  Step 

xxin  (points) 

Real 

Cross-Shore  Distances 

zzin  (points) 

Real 

Corresponding  Depths 

Output  Variables: 

nnn 

Integer 

Number  of  Points  in  Input  Depth  Array 

xdeltgr 

Real 

Self  Adjusting  Cross-Shore  Grid  Step 

xxl  (points) 

Real 

Adjusted  Cross-Shore  Distances  from 
Depth  Profile 

xxin  (points) 

Real 

Adjusted  Cross-Shore  Distances 

zzin  (points) 

Real 

Corresponding  Depths 

Local  Variables: 

None. 

Subroutines  Called  from  C  REGRID  (  ): 

C_FINE 

C_REGRID  ( )  Called  from  Subroutines: 

C_IN_DEP 
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Figure  10.  Subroutine  C  REGRID  Flowchart 


START  ) 

l 

/  Number  of  \ 
Points  in  Input 
Depth  Array 
_ Loop _ 


I 
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5.10  Subroutine  C  UN 


Subroutine  Call: 

C  UN  (dcall,  dcal2,  ndepth,  xxin,  zzin,  xdeltgr,  dstart) 

Summary: 

Subroutine  C  UN  converts  measurement  units  of  input  cross-shore  distances,  depth  arrays, 
starting  depth  and  the  grid  step  size  (Ax)  to  meters  for  internal  calculations. 


Real  Conversion  Factor  for  Cross-Shore  Distances 

Real  Conversion  Factor  for  Water  Depths 

Real  Input  Starting  Depth 

Integer  Number  of  Points  in  Input  Depth  Profile 

Real  Self  Adjusting  Cross-Shore  Grid  Step 

Real  Cross-Shore  Distances 

Real  Corresponding  Depths 

Output  Variables: 

dstart 
xdeltgr 
xxin  (points) 
zzin  (points) 

Local  Variables: 

I  Integer  Loop  Counter 

Subroutines  Called  from  C_UN  ( ):  None. 

C  UN  ( )  Called  from  Subroutines: 

C  IN  DEP 


Real  Input  Starting  Depth 

Real  Self  Adjusting  Cross-Shore  Grid  Step 

Real  Cross-Shore  Distances 

Real  Corresponding  Depths 


Input  Variables: 

dcall 

dcal2 

dstart 

ndepth 

xdelt_gr 

xxin  (points) 

zzin  (points) 
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Figure  11.  Subroutine  C  UN  Flowchart 
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5.11  Subroutine  CALC  HB3 


Subroutine  Call: 

CALC  HB3  (dp,  hrms,  p_flag,  hb3) 

Summary: 

Subroutine  CALC  HB3  integrates  the  wave  height  distribution  for  a  given  root  mean  square 
wave  height  and  calculates  a  term  inherent  to  the  roller  dissipation  function. 

Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

hrms 

Real 

Root  Mean  Square  Wave  Height  Calculation 

pflag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

Output  Variables: 

' 

hb3 

Real 

Weighting  Function  for  Dissipation  Term 

Local  Variables: 

hhigh 

Real 

Maximum  Wave  Height 

hlow 

Real 

Minimum  Wave  Height 

integrat 

Real 

Wave  Height  Distribution  Calculated  for  a 
Single  Wave  at  a  Specific  Location 

Functions  Called  from  CALC  HB3  (  ): 

INTEGRAT 


CALCJHB3  (  )  Called  from  Subroutines: 

CALCROLL 

GETDISS 
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Figure  12.  Subroutine  CALC_HB3  Flowchart 
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5.12  Subroutine  CALCROLL 


Subroutine  Call: 

CALCROLL  (roller,  hrms,  dp,  fqz,  theta,  xk,  e_roller) 

Summary: 

Subroutine  CALCROLL  calculates  roller  energy  at  a  point  in  the  surf  zone  based  on  water 
depth  and  Wave  Height  (hrms)  at  that  location. 


Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

fqz 

Real 

Zero  Crossing  Frequency 

hrms 

Real 

Root  Mean  Square  Wave  Height 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

theta 

Real 

Wave  Angle,  Representative  of  Radiation 
Stress  Angle 

xk 

Real 

Wave  Number 

Output  Variables: 

eroller 

Real 

Roller  Contribution  to  Energy  Equation 

Local  Variables: 

c 

Real 

Wave  Celerity 

er 

Real 

Temporary  Roller  Variable 

hb3 

Real 

Weighting  Function  for  Dissipation  Term 

pflag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

z 

Real 

Roller  Energy  Multiplier 

Subroutines  Called  from  CALCROLL(  ): 

CALCHB3 

CALCROLL( )  Called  from  Subroutines: 
GETRHS 
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CALCROLL( )  Called  from  Functions: 

F3 

Figure  13.  Subroutine  CALCROLL  Flowchart 


J 


RETURN  ◄  — 
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5.13  Subroutine  CALCSURF 


Subroutine  Call: 

CALCSURF  ( roller,  linstress,  ehsig,  wspd,  wdir,  tide,  ydepth,  nnn,  dxyl,  xxl,  ifreq,  freql, 
freq2,  freq,  idirec,  xfrom,  esowm,  dstart,  igamma.  ydetail,  iyear,  imonth,  iday,  ihour,  imin,  xdelt, 
xdelt_gr,  self_st,  surf,  pet,  alfa,  bravo,  chrlie,  echo,  foxtrt,  golfl,  golf2,  ihtll,  ihtl2,  jgamma) 

Summary: 

Subroutine  CALCSURF  acts  as  the  primary  driver  for  the  various  subroutines,  which 
calculate  wave  parameters  and  the  longshore  current  across  the  surf  zone. 


Input  Variables: 

dstart 

Real 

dxyl  (points) 

Real 

ehsig 

Real 

esowm  (dirNum,freqNum) 

Real 

freq  (freqNum) 

Real 

freql  (freqNum) 

Real 

freq2  (freqNum) 

Real 

iday 

Integer 

idirec 

Integer 

ifreq 

Integer 

igamma 

Integer 

ihour 

Integer 

imin 

Integer 

imonth 

Integer 

iyear 

Integer 

lin_stress 

Logical 

nnn 

Integer 

roller 

Logical 

selfst 

Char*  1 

tide 

Real 

wdir 

Real 

wspd 

Real 

xdelt 

Real 

xdelt_gr 

Real 

xfrom  (dirNum) 

Real 

Input  Starting  Depth 
Corresponding  Depths  without  Tide 
Significant  Wave  Height  from 
Directional  Spectrum 
Directional  Wave  Spectrum 
Input  Wave  Spectrum  Center  Frequency 
Beginning  Frequency  Bin  Value 
Ending  Frequency  Bin  Value 
Input  Day 

Number  of  Direction  Bins  in  Input  Spectrum 
Number  of  Frequencies  in  Input  Spectrum 
Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 
Input  Hour 
Input  Minute 
Input  Month 
Input  Year 

Longshore  Current  Solution  (True  or  False) 
Number  of  Points  in  Input  Depth  Array 
Roller  Option  Flag  (True  or  False) 

Self  Start  Flag  (Yes  or  No) 

Input  Tide  Level 

Input  Wind  Direction,  Compass  Heading  Wind 

is  Blowing  From 

Input  Wind  Speed 

Surf  Zone  Output  Interval 

Self-Adjusting  Cross-Shore  Grid  Step 

Direction  Array,  Direction  Wave  Energy 
Comes  From 
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xx  1  (points) 

Real 

ydepth 

Char*  1 

ydetail 

Char*  1 

Output  Variables: 

alfa 

Real 

bravo 

Real 

chrlie 

Real 

echo 

Real 

foxtrt 

Real 

golfl 

Real 

golf2 

Real 

ihtll 

Real 

ihtl2 

Real 

j  gamma 

Integer 

pet (4) 

Real 

surf 

Logical 

Local  Variables: 

along  (points) 

Real 

b 

Real 

bl  (points) 

Real 

blong  (points) 

Real 

c 

Real 

cl 

Real 

c2 

Real 

c3 

Real 

c4 

Real 

cf 

Real 

Cg 

Real 

clong  (points) 

Real 

convg 

Logical 

df 

Real 

Adjusted  Cross-Shore  Distances  from 
Depth  Profile 

Input  Depth  Profile  Used?  (Yes  or  No) 
Detailed  Output?  (Yes  or  No) 


Significant  Breaker  Height 
Maximum  Breaker  Height 
Dominant  Breaker  Period 
Breaker  Angle 

Longshore  Current  Speed  and  Direction 

Number  of  Surf  Lines 

Surf  Zone  Width 

Wind  Speed 

Wind  Direction 

Temporary  Value  Set  to  Beach  Onentation 

Percent  of  Different  Breaker  Types: 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

Flag  for  Low/No  Surf  Conditions 
(True  or  False) 


Horizontal  Mixing  Parameter  from 
Thornton  &  Whittord 

Empirical  Factor  in  Thornton  &  Guza  Wave 
Breaking  Model  (=  1 .00) 

Bottom  Slope 

Bottom  Friction  for  Deep  &  Shallow  Water 

Wave  Celerity  at  Input  Starting  Depth 

Eddy  Viscosity  Coefficient 

Bottom  Friction  Coefficient 

Radiation  Stress  Coefficient  -  Multiple  for 

Longshore  Current  Model 

Longshore  Wind  Stress  Coefficient  -  Multiple 

for  Longshore  Current  Model 

Coefficient  of  Bottom  Friction 

Wave  Group  Velocity 

Wind  Stress  Contribution  to 

Longshore  Current 

Energy  Equation  Convergence  Flag 

Difference  Between  Adjacent  Frequency  Bins 
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distmax 

Real 

dp 

Real 

dth 

Real 

dwsstop 

Integer 

dxy  (points) 

Real 

eb_last 

Real 

ebtemp  (points) 

Real 

filespc 

Char*40 

fqd 

Real 

fqz 

Real 

ftsq2msq 

Real 

hi  max 

Real 

h2max 

Real 

hrms 

Real 

htemp  (points) 

Real 

iimax 

Integer 

irealf 

Integer 

j 

Real 

jji 

Integer 

j_ii2 

Integer 

k 

Real 

per 

Real 

print_spc 

Integer 

ptemp  (points) 

Real 

rk  (points, 4) 

Real 

stringout 

Character 

stringsub 

Character 

suml 

Real 

temp 

Real 

theta 

Real 

thetal 

Real 

Farthest  Offshore  Distance 
Offshore  Water  Depth 

Difference  Between  Adjacent  Directional  Bins 
Flag  for  Shallow  Water  Directional  Wave 
Spectrum  Print  Control 
Pre-Tidal  Depth  with  Tide 
Roller  Dissipation  Term  at  Farthest 
Point  Offshore 

Temporary  Roller  Dissipation  Term 
Across  Transect 

File  Name  of  Shallow  Water  Directional  Wave 
Spectrum 

Peak  Frequency  at  the  Center  of  the  Frequency 
Band 

Zero  Crossing  Frequency 

Conversion  Factor  from  Feet  Squared  to  Meters 

Squared 

Largest  Significant  Wave  Height  in  the 
Surf  Zone 

Largest  Maximum  Wave  Height  in  the 
Surf  Zone 

Root  Mean  Square  Wave  Height 
Temporary  Variable  for  Significant  Wave 
Height  Values 

Number  of  Calculation  Locations 
Cutoff  Index  for  Printing  Shallow  Water 
Directional  Wave  Spectrum 
Temporary  Variable  for  Cross-Shore  Values 
Index  where  Wave  Probabilities  come 
Above  Threshold 

Longshore  Current  Loop  Variable  for  Outer 
Edge  of  Surf  Zone 
Temporary  Variable  for  Significant 
Wave  Height 

Peak  Period  of  Directional  Wave  Spectrum 
Flag  to  Print  Shallow  Water  Wave  Spectrum 
Percentage  of  Breaking  Waves  and 
Breaker  Types 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 

Shallow  Water  Wave  Spectrum  Output  String 

Temporary  String  Variable 

Sum  of  Wave  Length  in  the  Surf  Zone 

Temporary  Variable 

Wave  Angle 

Wave  Angle  at  Input  Starting  Depth 
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theta2 

Real 

Wave  Angle  at  Input  Starting  Depth 

v  (points) 

Real 

Longshore  Current 

vmax 

Real 

Maximum  Positive  Longshore  Current 

vmin 

Real 

Maximum  Negative  Longshore  Current 

vwind 

Real 

Group  Wind  Velocity 

wdspd 

Real 

Wind  Speed  Conversion 

Knots  to  CM/S  =  5 1 .44 

widii 

Integer 

Array  Location  for  Surf  Zone  Width 

width 

Real 

Surf  Zone  Width 

xk 

Real 

Wave  Number 

xktemp  (points) 

Real 

Temporary  Variable  for  Wave  Number 

xshift 

Real 

Horizontal  Cross-Shore  Location 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

Subroutines  Called  from  CALCSURF  ( ): 

BDETAIL 

GTSIGH 

INITLIZE 

KLONG 

MAIN_WAV 

NONLIN 

PRT_OUTl 

PRTOUT2 

RADST1 

SCOEFF 

SNOSURF 

STEDE 

SHORTOUT 

ZONE1 

CALCSURF  ( )  Called  from  Subroutines: 

SURF 
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5.14  Subroutine  CONANG 
Subroutine  Call: 

CON_ANG  (t,  h,l,  dp,  q,  theta2,  u,  v,  convg,  kount) 


Summary: 

Subroutine  CON  ANG  calculates  the 


momentum  in  the  longshore  direction. 

Input  Variables: 


dp 

Real 

h 

Real 

1 

Real 

q 

Real 

t 

Real 

theta2 

Real 

u 

Real 

Output  Variables: 

convg 

Logical 

kount 

Integer 

V 

Real 

Local  Variables: 

fl 

Real 

£2 

Real 

numit 

Integer 

tol 

Real 

v  new 

Real 

longshore  current  velocity  based  on  the  flux  of 


Offshore  Water  Depth 
Wave  Height 
Wave  Length 

Longshore  Momentum  Flux 
Wave  Period 
Rotated  Wind  Direction 
Cross-Shore  Current  Velocity 


Convergence  Flag  (True  or  False) 
Counter 

Longshore  Current  Velocity 


Wave  Height  Distribution 
Weighting  Function 
Wave  Height  Distribution 
Weighting  Function 
Number  Limitation  -  Set  to  1000 
Tolerance  Check  -  Set  to  1.0E-4 
Temporary  Longshore  Current  Velocity 


Subroutines  Called  from  CON  ANG  ( ): 

GETFCN 
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CONANG  ( )  Called  from  Subroutines: 

NONLIN2 

Figure  15.  Subroutine  CON  ANG  Flowchart 


5.15  Subroutine  DEPDRVR 


Subroutine  Call: 

DEPDRVR  (depname,  dstart,  xdelt,  ydepth,  slope,  gamma2,  nnn,  xxl,  dxyl,  igamma, 
xdelt_gr) 

Summary: 

Subroutine  DEPDRVR  is  the  driver  routine  for  reconstructing  the  depth  arrays  in  an 
optimized  step  size. 


Input  Variables: 

depname 

Char*40 

dstart 

Real 

gamma2 

Real 

slope 

Real 

xdelt 

Real 

ydepth 

Char*  1 

Output  Variables: 

dxyl  (points) 

Real 

igamma 

Integer 

nnn 

Integer 

xdelt_gr 

Real 

xxl  (points) 

Real 

Local  Variables: 

None. 

Subroutines  Called  from  DEPDRVR  ( ): 

EQUILPRF 

CGAMMA 

CJNJDEP 


DEPDRVR  ( )Called  from  Subroutines: 

SURF 

Figure  16.  Subroutine  DEPDRVR  Flowchart 


Depth  Profile  File  Name 

Input  Starting  Depth 

Beach  Orientation  Compass  Heading 

Directly  Toward  Beach 

Bottom  Slope 

Surf  Zone  Output  Interval 

Usage  of  Input  Depth  Profile  (Yes  or  No) 


Corresponding  Depths  without  Tide 
Beach  Orientation  Rotated  90  Degrees  from  the 
Original  Heading  Toward  the  Beach 
Number  of  Points  in  the  Input  Depth  Array 
Self-Adjusting  Cross-Shore  Grid  Step 
Adjusted  Cross-Shore  Distances  from  the  Depth 
Profile 
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5.16  Subroutine  EQUILPRF 
Subroutine  Call: 

EQUILPRF  (rtype,  dpthoff,  xgrd,  numstep,  xxl,  dxyl) 

Summary: 


Subroutine  EQUILPRF  constructs  a  depth  profile  for  surf  calculations.  This  equilibrium 
profile  is  based  on  the  equation:  y=AxA(2/3),  where  A  is  a  coefficient  related  to  sediment  grain  size 
or  frictional  dissipation.  This  equation  was  developed  by  Dean  (1977)  from  a  study  of  more  than  200 
beach  profiles.  The  "A"  coefficient  in  the  equilibrium  equation  has  units  of  meters,  calculations  in 
feet  require  different  values  or  conversion  to  feet  after  initial  calculations.  Sediment/grain  types  are 
denoted  by  the  variable  "rtype"  which  is  the  index  for  a  value  in  the  array  of  coefficients  defining  the 
following  grain  sizes: 

1  =  boulders 

2  =  cobble 

3  =  pebbles 

4  =  granules 

5  =  very  coarse  sand 

6  =  coarse  sand 

7  =  medium  sand 

8  =  fine  sand 

9  =  very  fine  sand 

10  =  silt 


Input  Variables: 

dpthoff 

Real 

Input  Starting  Depth 

numstep 

Integer 

Number  of  Points  in  the  Input  Depth  Array 

rtype 

Real 

Sediment/  Grain  Type 

xgrd 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

Output  Variables: 

dxyl  (points) 

Real 

Corresponding  Depths  with  No  Tide 

xxl  (points) 

Local  Variables: 

Real 

Cross-Shore  Distances 
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a(10) 

Real 

ause 

Real 

call 

Real 

distance 

Logical 

diston 

Real 

dpthon 

Real 

I 

Integer 

X 

Real 

xone 

Real 

z 

Real 

Subroutines  Called  from  EQUILPRF  (  ): 
EQUILPRF  ( )  Called  from  Subroutines: 

DEPDRVR 


Array  of  Sediment  Coefficients 

Actual  Sediment  Type  Coefficient  for  Profile 

Conversion  Factor  (Meters) 

Flag  for  Equilibrium  Depth  Bottom 
Highest  Onshore  Distance 
Highest  Onshore  Depth 
Loop  Counter 
Temporary  Variable 
Farthest  Point  Offshore 
Temporary  Variable 

ABORT 
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Figure  17.  Subroutine  EQUILPRF  Flowchart 
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5.17  Subroutine  GENRLSPL 


Subroutine  Call: 

GENRLSPL  (xin,  yin,  inlen,  xout,  outlen,  yout) 

Summary: 

Subroutine  GENRLSPL  is  the  driver  routine  to  interpolate  an  array  of  x  and  y  values  to  a  new 


set  of  x  values  using  a  cubic  spline  polynomial. 

Input  Variables: 


inlen  Integer 

outlen  Integer 

xin  (dirNum)  Real 

xout  (dirNum)  Real 

yin  (dirNum)  Real 

Output  Variables: 

yout  (dirNum)  Real 

Local  Variables: 

coef  (4,dirNum)  Real 

cubpoly  Real 

I  Integer 


Subroutines  Called  from  GENRLSPL  ( ): 

SPLINE 

Functions  Called  from  GENRLSPL  (  ): 

CUBPOLY 

GENRLSPL  (  )  Called  from  Subroutines: 

READRFRC 


Number  of  input  Coordinates 
Number  of  Coordinates  to  Interpolate 
X-Coordinates  of  known  Values 
Interpolated  X-Coordinates 
Y-Coordinates  of  known  Values 


Interpolated  Y-Coordinates 


Temporary  Array  of  Interpolated  Coefficients 
Value  at  the  Interpolated  Coordinate 
Loop  Counter 
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Figure  18.  Subroutine  GENRLSPL  Flowchart 


I 
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5.18  Subroutine  GENSPEC 


Subroutine  Call: 

GENSPEC  (hsea,  psea,  dsea,  hswell,  pswell,  dswell,  ifreq,  idirec,  freq,  freql,  freq2,  xfrom, 
esowm,  period,  ehsig,  dangle) 

Summary: 

Subroutine  GENSPEC  initializes  matrices  for  the  creation  of  an  internally  generated 


directional  wave  spectrum. 

This  wave  spectrum  has  50  frequencies  and  36  directions. 

Input  Variables: 

dsea 

Real 

Input  Direction  for  Sea  Contribution 

dswell 

Real 

Input  Swell  Direction  for  Internally 
Generated  Spectrum 

hsea 

Real 

Input  Significant  Wave  Height  for  Sea 
Contribution  to  Pierson  Moskowitz  Equation 

hswell 

Real 

Input  Significant  Wave  Height  for  Internally 
Generated  Spectrum 

psea 

Real 

Input  Wave  Period  for  Sea 

Contribution  to  Pierson  Moskowitz  Equation 

pswell 

Real 

Input  Swell  Period  for  Internally 

Generated  Spectrum 

Output  Variables: 

dangle 

Real 

Angle  Between  Directional  Bins 

ehsig 

Real 

Significant  Wave  Height  from 

Directional  Spectrum 

esowm  (dirNum,ffeqNum) 

Real 

Directional  Spectrum 

freq  (freqNum) 

Real 

Input  Wave  Spectrum  Center  Frequencies 

freql  (freqNum) 

Real 

Beginning  Frequency  Bin  Values 

freq2  (freqNum) 

Real 

Ending  Frequency  Bin  Values 

idirec 

Integer 

Number  of  Direction  Bins  in  the 

Input  Spectrum 

ifreq 

Integer 

Number  of  Frequencies  in  the 

Input  Spectrum 

period  (freqNum) 

Real 

Period  Array  (1/Frequency) 

xfrom  (dirNum) 

Real 

Direction  Array,  Direction  Wave  Energy 
Comes  From 
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Local  Variables: 


df 

Real 

Difference  between  Frequency  Bins 

idir 

Integer 

Direction  Loop  Counter 

ifrq 

Integer 

Frequency  Loop  Counter 

Subroutines  Called  from  GENSPEC  ( ): 
WAVEFIT 

GENSPEC  ( )  Called  from  Subroutines: 

SURF 

Figure  19.  Subroutine  GENSPEC  Flowchart 


[ START 

I 

Set  Frequency 
&  Direction 
Limits 


Initialize 
Frequency  Bins 


Initialize 
Period  .Array 

~r 


End  of  Loop 


l 


Call 

WAVEFIT 


I 

(  RETURN  ] 
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5.19  Subroutine  GETBRK 
Subroutine  Call: 

GET  BRK  ( ii,  dxy,  xdelt_gr,  hrms,  per,  xoff,  rk,  bl,  brklO,  distmax,  p) 

Summary: 

Subroutine  GET  BRK  calculates  percentage  of  breakers  and  percent  breaker  type  given  at 
each  point  along  the  transect:  p  (1)  =  Spilling,  p  (2)  =  Plunging,  p  (3)  =  Surging,  p  (4)  =  100*Sum. 


Input  Variables: 


bl  (points) 

Real 

Bottom  Slope 

brklO 

Logical 

Flag  for  First  Location  where  10%  of  the 
Waves  are  Breaking  (True  or  False) 

distmax 

Real 

Farthest  Offshore  Distance 

dxy  (points) 

Real 

Adjusted  Depths  with  Tide 

hrms 

Real 

Root  Mean  Square  Wave  Height 

ii 

Integer 

Index  where  Wave  Probabilities  Exceed 
Threshold 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

rk  (points,4) 

Real 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 

xdelt_gr 

Real 

Self-Adjusting  Cross-Shore  Step 

xoff 

Real 

Distance  Offshore 

Output  Variables: 

bl  (points) 

Real 

Bottom  Slope 

brklO 

Logical 

Flag  for  First  Location  where  1 0%  of  the  Waves 
are  Breaking  (True  or  False) 

distmax 

Real 

Farthest  Offshore  Distance 

P(4) 

Real 

Temporary  Array  for  Breaker 

Percentage  Totals 

rk  (points, 4) 

Real 

Percent  Breaker  of  Each  Type 

Local  Variables: 

beta 

Real 

Temporary  Variable  for  Bottom  Slope 
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Subroutines  Called  from  GETBRK  (  ): 

PERCENT 


GET_BRK  ( )  Called  from  Subroutines: 

MAINWAV 

Figure  20.  Subroutine  GET  BRK  Flowchart 

[  START  ] 


I 


I 


Set  Breaking  Wave 
Percentage  Flags 


I 

[  RETURN 
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5.20  Subroutine  GETDISS 


Subroutine  Call: 

GET  DISS  (roller,  b,  fqz,  dp,  hrms,  pjlag,  diss) 

Summary: 

Subroutine  GET  DISS  returns  the  wave  dissipation  factor.  This  term  is  based  on  a  Bore 
dissipation  Model  and  can  include  roller  dissipation  if  selected.  The  dissipation  term  is  included  in 


the  wave  energy  balance  equation.  The  wave  dissipation  is  given  by: 

Where  <p  is  density,  g  is  gravity,  f  is  bottom  friction,  h  is  the  water  depth,  M  is  a  weighting  function 
based  on  hrms,  and  B  is  an  empirical  factor. 


Input  Variables: 


b 

Real 

Empirical  Factor  in  Thornton  &  Guza  Wave 
Breaking  Model  =  1 .00 

dp 

Real 

Offshore  Water  Depth 

fqz 

Real 

Zero  Crossing  Frequency 

hrms 

Real 

Root  Mean  Square  Wave  Height 

p _ flag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

Output  Variables: 

diss 

Real 

Bore  or  Roller  Dissipation  Function 
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Local  Variables: 


Real 

Real 


Weighting  Function  for  Dissipation  Term 
Dissipation  Function 


hb3 

z 

Subroutines  Called  from  GET_DISS  ( ): 

CALCHB3 

GETDISS  ( )  Called  from  Subroutines: 

GETRHS 

Figure  21.  Subroutine  GET  DISS  Flowchart 
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5.21  Subroutine  GETFCN 
Subroutine  Call: 

GET  FCN  (t,  h,  1,  dp,  v,  u,  theta2,  fl,  f2) 

Summary: 

Subroutine  GET  FCN  performs  a  call  to  Function  FCN1  to  evaluate  the  nonlinear  bottom 
stress  at  a  specific  time  interval.  This  function  integrates  the  bottom  stress  over  time  utilizing  a 
trapezoidal  integration  method  to  provide  the  time-averaged  bed  stress  at  a  certain  location. 


Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

h 

Real 

Wave  Height 

1 

Real 

Wave  Length 

t 

Real 

Wave  Period 

theta2 

Real 

Wave  Angle 

u 

Real 

Cross-Shore  Current  Velocity 

V 

Real 

Longshore  Current  Velocity 

Output  Variables: 

fl 

Real 

Wave  Height  Distribution  Weighting  Function 

f2 

Real 

Wave  Height  Distribution  Weighting  Function 

Local  Variables: 

delt 

Real 

Time  Step  Interval 

f_xo_l 

Real 

Integral  Evaluated  at  the  Lower 

fxnl 

Real 

Limit  of  Integration 

Integral  Evaluated  at  Upper  Limit 

fcnl 

Real 

of  Integration 

Nonlinear  Bottom  Friction  at  a  Specific  Time 

i 

Integer 

Loop  Increment 

numit 

Integer 

Set  Equal  to  100 

suml 

Real 

Local  Integration  Variable 

xi 

Real 

Integration  Step  Location 

xn 

Real 

Upper  Limit  of  Integration 

xo 

Real 

Lower  Limit  of  Integration 
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Functions  Called  from  GETFCN  (  ): 

FCN1 

GETJFCN  ( )  Called  from  Subroutines: 

CONANG 


Figure  22.  Subroutine  GET  FCN  Flowchart 
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5.22  Subroutine  GETM 


Subroutine  Call: 

GET_M  (dp,  hrms,  m) 

Summary: 

Subroutine  GETM  calculates  the  weighting  function  multiplier. 

Input  Variables: 

dp  Real  Offshore  Water  Depth 

hrms  Real  Root  Mean  Square  Wave  Height 

Output  Variables: 

m  Real  Multiplier 

Local  Variables:  None. 

Subroutines  Called  from  GET_M  (  ):  None. 

GET  M  ( )  Called  from  Subroutines: 

WEIGHTFN 

Figure  23.  Subroutine  GET  M  Flowchart 


"  START  ] 

\ 

Calculate  Multiplier 
Based  Upon 
Weighting  Function 

i 

RETURN I 
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5.23  Subroutine  GET  P 


Subroutine  Call: 

GET_P  (frac,  p) 

Summary: 

Subroutine  GETJP  calculates  the  percentage  of  each  breaker  type  and  fills  the  corresponding 

array  elements. 

Input  Variables: 

frac  (3)  Real  Temporary  Array  for  Breaker 

Percentage  Totals 

Output  Variables: 

P  (4)  Real  Percent  of  Different  Breaker  Types 

P  (1)  =  Spilling 
p  (2)  =  Plunging 
p  (3)  =  Surging 
p  (4)  =  Total 

Local  Variables: 

sum  Real  Temporary  Variable  for  Total  of 

Percentage  Breakers 

Subroutines  Called  from  GET  P  ( ):  None. 

GET_P  ( )  Called  from  Subroutines: 

PERCENT 
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Figure  24.  Subroutine  GET  P  Flowchart 
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5.24  Subroutine  GETRHS 
Subroutine  Call: 

GET  RHS  (roller,  hrms,  theta,  Cg,  dp,  xk,  b,  fqz,  xdelt_gr,  rhs,  diss) 

Summary: 

Subroutine  GET  RHS  calculates  the  right  hand  side  of  the  wave  energy  equation. 

Input  Variables: 


b 

Real 

Empirical  Factor  in  Breaking  Model  =  1 .0 

Cg 

Real 

Wave  Group  Velocity 

dp 

Real 

Offshore  Water  Depth 

fqz 

Real 

Zero  Crossing  Frequency 

hrms 

Real 

Root  Mean  Square  Wave  Height 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

theta 

Real 

Wave  Angle,  Representative  of  Radiation 
Stress  Angle 

xdelt_gr 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

xk 

Real 

Wave  Number 

Output  Variables: 

diss 

Real 

Bore  or  Roller  Dissipation  Function 

rhs 

Real 

Right  Hand  Side  of  Energy  Equation 

Local  Variables: 

e_roller 

Real 

Roller  Contribution  to  the  Energy  Equation 

e_wave 

Real 

Wave  Contribution  to  the  Energy  Equation 

pflag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

Subroutines  Called  from  GET  RHS  ( ): 


CALCROLL 

GETDISS 

GET_WAVE 

GET_RHS  (  )  Called  from  Subroutines: 

MAIN_WAV 


Figure  25.  Subroutine  GET  RHS  Flowchart 
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5.25  Subroutine  GET  UM 


Subroutine  Call: 

GET_UM  (h,  t,  1,  dp.  urn) 

Summary: 

Subroutine  GETUM  uses  linear  wave  theory  to  calculate  the  wave-induced  orbital  velocity. 
The  wave  induced  orbital  velocity  is  calculated  where: 

gHT 

u - - ~2tc 

2  L  cosh  ( — )  h 
L 


g  is  gravity,  H  is  wave  height.  T  is  wave  period,  L  is  wave  length,  and  h  is  water  depth. 


Input  Variables: 

dp 

Real 

Offshore  Water  Depth 

h 

Real 

Wave  Height 

1 

Real 

Wave  Length 

t 

Real 

Wave  Period 

Output  Variables: 

um 

Real 

Wave  Induced  Orbital  Velocity 

Local  Variables: 

None. 

Subroutines  Called  from  GET  UM  (  ): 

None. 

GET_UM  (  )  Called  from  Functions: 
FCN1 
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Figure  26.  Subroutine  GETJLJM  Flowchart 


START 

T 


Calculate  Wave 
Orbital  Velocity 


_L 

RETURN 
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5.26  Subroutine  GETWAVE 
Subroutine  Call: 

GET_WAVE  (hrms,  theta,  Cg,  e  wave) 

Summary: 

Subroutine  GET  WAVE  calculates  wave  energy  flux  using  linear  wave  theory.  The  wave 
energy  flux  is: 


E 


<Pg  //: 
8 


Cg  cos  6 


where  q>  is  water  density,  g  is  gravity,  H  is  wave  height,  Cg  if  group  velocity,  and  0  is  the  wave  angle. 

Input  Variables: 

Cg  Real  Wave  Group  Velocity 

hrms  Real  Root  Mean  Square  Wave  Height 

theta  Real  Wave  Angle 

Output  Variables: 


e  wave  Real  Energy  Flux 

Local  Variables: 


ew  Real  Wave  Energy 

Subroutines  Called  from  GET  WAVE  ( ):  None. 


GET  WAVE  ( )  Called  from  Subroutines: 


F3 

GETRHS 
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Figure  27.  Subroutine  GET  WAVE  Flowchart 
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5.27  Subroutine  GRIDOUT 


Subroutine  Call: 


GRIDOUT  (ii,  xoffl,  xtemp,  dxy,  htemp,  ptemp,  xktemp,  v,  dpi,  houtl,  hmax,  pbreak, 
wlen,  vlngl) 

Summary: 

Subroutine  GRIDOUT  linearly  interpolates  parameters  for  final  output  using  the  user  defined 
cross-shore  step  width. 

Input  Variables: 


dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

ii 

Integer 

Index  where  Wave  Probabilities  Exceed 
Threshold 

ptemp  (points) 

Real 

Percentage  of  Breaking  Waves  and 

Breaker  Types 

xktemp 

Real 

Temporary  Variable  for  Wave  Number 

xoffl 

Real 

Distance  Offshore 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

v  (points) 

Real 

Longshore  Current 

Output  Variables: 

dpi 

Real 

Offshore  Depth 

hmax 

Real 

Maximum  Wave  Height  /  10.00 

houtl 

Real 

Significant  Wave  Height 

pbreak 

Real 

Percentage  Breaking  Waves 

vlngl 

Real 

Longshore  Current  Velocity 

wlen 

Real 

Wave  Length 

Local  Variables: 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 
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Subroutines  Called  from  GRIDOUT  ( ): 

GT_P 

LINl 

LIN_2 

LIN3 

GRIDOUT  ( )  Called  from  Subroutines: 

PRTOUT1 

Figure  28.  Subroutine  GRIDOUT  Flowchart 

[  START  ) 


I 
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5.28  Subroutine  GT  P 


Subroutine  Call: 


GT_  P  (ii,  hmisl,  dpi,  xktenip,  houtl,  hmax,  wlen) 

Summary: 


Subroutine  GT_P  initializes  matrices  for  the  creation  of  an  internally  generated  directional 
wave  spectrum.  This  wave  spectrum  has  50  frequencies  and  36  directions. 


Input  Variables: 

dpi  Real 

ii  Integer 

hrmsl  Real 

xktemp  (points)  Real 

Output  Variables: 

hmax  Real 

houtl  Real 

wlen  Real 

Local  Variables: 

hdep  Real 


Subroutines  Called  from  GT  P  ( ): 
GT_P  (  )  Called  from  Subroutines: 


Offshore  Depth 

Index  where  Wave  Probabilities 
Exceed  Threshold 
Root  Mean  Square  Wave  Height 
Temporaiy  Variable  for  Wave  Number 


Maximum  Wave  Height  /  JO.OO 
Significant  Wave  Height 
Wave  Length 


Breaking  Wave  Height  at  Specific  Depth 
None. 


GRIDOUT 

PRTOUT1 

PRT_OUT2 
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Figure  29.  Subroutine  GT_P  Flowchart 


START 

T 


Calculate  Wave  Heights 
and  Breaker  Height 


_L 

RETURN 
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5.29  Subroutine  GT  SIG  H 


Subroutine  Call: 

GT_SIG_H  (ifreq,  idirec,  esowm,  ehsig) 

Summary: 

Subroutine  GTSIGH  sums  the  energy  present  in  the  directional  wave  spectrum  and 


jJZetfO) 


calculates  the  significant  wave  height.  The  significant  wave  height  is  defined  as: 
Where,  e  is  the  directional  wave  spectrum. 


Input  Variables: 

esowm  (dirNum,freqNum)  Real 
idirec  Integer 

ifreq  Integer 

Output  Variables: 

ehsig  Real 


Local  Variables: 

idir  Integer 

ifrq  Integer 

suml  Real 

sum2  Real 

Subroutines  Called  from  GT  SIG_H  ( ): 

GT  SIG  H  ( )  Called  from  Subroutines: 

CALCSURF 

READSPEC 

WAVEFIT 


Directional  Wave  Spectrum 

Number  of  Direction  Bins  in  Input  Spectrum 

Number  of  Frequencies  in  Input  Spectrum 


Significant  Wave  Height  from 
Directional  Spectrum 


Direction  Loop  Counter 
Frequency  Loop  Counter 
Summing  Variable  for  Wave  Height 
Summing  Variable  for  Wave  Height 

None. 
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Figure  30.  Subroutine  GT_SIG_H  Flowchart 


I 
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Significant  Wave 
Height 

\ 

RETURN  ] 
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5.30  Subroutine  INITLIZE 


Subroutine  Call: 

INITLIZE  (dp,  fqd,  Cg,  xk,  c) 

Summary: 

Subroutine  INITLIZE  calculates  wave  parameters  at  the  farthest  offshore  point.  Wave  celerity 


o2  =  gk  tanh (k hj 

(velocity)  is  calculated  from  the  dispersion  relation  given  by: 

where,  a  is  the  angular  frequency  of  the  wave  (2xT),  g  is  gravity,  k  is  wave  number,  and  h  is  the  local 


Ce 


0.501  + 


2kh 
sinh  k  h 


water  depth.  Wave  group  velocity  is  calculated  from  the  linear  wave  theory  relation  given  by: 
where,  C  is  the  wave  celerity. 

Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

fqd 

Real 

Peak  Frequency 

Output  Variables: 

c 

Real 

Wave  Celerity  at  Input  Depth  &  Frequency 

Cg 

Real 

Group  Velocity  at  Input  Depth  &  Frequency 

xk 

Real 

Wave  Number  at  Input  Depth  &  Frequency 
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Local  Variables: 


xkd  Real  Deep  Water  Wave  Number 

Subroutines  Called  from  INITLIZE  ( ): 

WAVENUM 

INITLIZE  (  )  Called  from  Subroutines: 

CALCSURF 

Figure  31.  Subroutine  INITLIZE  Flowchart 


I 
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5.31  Subroutine  KLONG 


Subroutine  Call: 

KLONG  (j_ii,  xdelt_gr,  eb_last,  along,  blong,  clong,  c3,  iimax,  vwind,  v) 

Summary: 

Subroutine  KLONG  calculates  longshore  current  velocity  using  an  implicit  double  sweep 
method  (Tridiagonal  Method)  based  on  the  work  of  Kraus  and  Larson  (1991).  The  central  difference 


a,  Vj.,+  bi  Vr  a  Vi+i  =  n 


equation  is  of  the  form: 

where,  V  is  the  longshore  current  velocity.  The  coefficients  a,  b,  and  c  are  calculated  from  wave 
parameters. 

Input  Variables: 


along  (points) 

Real 

Horizontal  Mixing  Parameter 

blong  (points) 

Real 

Bottom  Friction 

c3 

Real 

Radiation  Stress  Factor  for  Longshore  Current 
Velocity 

clong  (points) 

Real 

Wind  Stress  Contribution  to 

Longshore  Current 

eblast 

Real 

Roller  Dissipation  Term  Farthest  Offshore 

iimax 

Integer 

Number  of  Calculation  Locations 

j_ii 

Integer 

Index  where  Wave  Probabilities 

Exceed  Threshold 

vwind 

Real 

Wind  Driven  Longshore  Current  Velocity 

xdeltgr 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

Output  Variables: 

v  (points) 

Real 

Longshore  Current  Velocity 
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Local  Variables: 


ah 

Real 

Temporary  Variable  Used  in  Longshore  Current 
Calculation 

bh 

Real 

Temporary  Variables 

ch 

Real 

Temporary  Variables 

dn 

Real 

Temporary  Variables 

ee  (points) 

Real 

Array  of  Longshore  Driving  Terms 

ff  (points) 

Real 

Array  of  Longshore  Bottom  Friction 

ieeff 

Integer 

Array  Index 

ii 

Integer 

Loop  Variable 

iuse 

Integer 

Array  Index  /  Loop  Variable 

xde!2 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

Subroutines  Called  from  KLONG  (  ): 

None. 

KLONG  (  )  Called  from  Subroutines: 

CALCSURF 
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Figure  32.  Subroutine  KLONG  Flowchart 
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5.32  Subroutines  LINl 
Subroutine  Call: 

LIN_1  (ii,  dx,  dy,  x,  y) 

Summary: 

Linear  interpolation  routine  used  to  scale  root  mean  square  wave  height  and  water  depth  to 
user-defined  grid  step  for  output  to  the  summary  text  file. 


Input  Variables: 


dx  (points) 

Real 

Input  X  Value 

dy  (points) 

Real 

Input  Y  Value 

ii 

Integer 

Index  where  Wave  Probabilities 
Exceed  Threshold 

X 

Real 

Offshore  Point 

Output  Variables: 

y 

Real 

Interpolated  Variable 

Local  Variables: 

bl 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yi 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_1  ( ):  None. 


LIN  l  ( )  Called  from  Subroutines: 

GRIDOUT 
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Figure  33.  Subroutine  LIN  l  Flowchart 
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5.33  Subroutines  LIN  2 
Subroutine  Call: 

LIN  2  (ii,  dx,  dy,  x,  y) 

Summary: 

Linear  interpolation  routine  used  to  scale  percent  breaking  waves  to  user-defined  grid  step  for 
output  to  the  summary  text  file. 


Input  Variables: 


dx  (points) 

Real 

Input  X  Value 

dy  (points) 

Real 

Input  Y  Value 

ii 

Integer 

Index  where  Wave  Probabilities 
Exceed  Threshold 

X 

Real 

Offshore  Point 

Output  Variables: 

y 

Real 

Interpolated  Variable 

Local  Variables: 

bl 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yi 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_2  ( ):  None. 


LIN  2  ( )  Called  from  Subroutines: 

GRIDOUT 
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Figure  34.  Subroutine  LIN  2  Flowchart 
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5.34  Subroutine  LIN  3 


Subroutine  Call: 

LIN_3  (ii,  dx,  dy,  x,  y) 

Summary: 

Linear  interpolation  routine  used  to  scale  longshore  current  velocity  distribution  to  user- 
defined  grid  step  for  output  to  the  summary  text  file. 


Input  Variables: 


dx  (points) 

Real 

Input  X  Value 

dy  (points) 

Real 

Input  Y  Value 

ii 

Integer 

Index  where  Wave  Probabilities 
Exceed  Threshold 

X 

Real 

Offshore  Point 

Output  Variables: 

y 

Real 

Interpolated  Variable 

Local  Variables: 

bl 

Real 

Intercept 

m 

Real 

Slope 

xl 

Real 

Cross-Shore  Value 

x2 

Real 

Previous  Cross-Shore  Value 

yl 

Real 

Height  Value 

y2 

Real 

Previous  Height  Value 

Subroutines  Called  from  LIN_3  (  ):  None. 


LIN  3  ( )  Called  from  Subroutines: 

GRIDOUT 
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Figure  35.  Subroutine  LIN  3  Flowchart 
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5.35  Subroutine  LONG1 


Subroutine  Call: 

LONG1  (ii,  cl,  c2,  c3,  c4,  dp,  ebn,  hrms,  xk,  along,  blong,  clong) 

Summary: 

Subroutine  LONG1  calculates  and  outputs  longshore  current  equation  coefficients. 

Input  Variables: 

cl 

c2 
c3 
c4 
dp 
ebn 
ii 

hrms 


Output  Variables: 


along  (points) 

Real 

Horizontal  Mixing  Parameter 

blong  (points) 

Real 

Bottom  Friction  Parameter 

clong  (points) 

Real 

Wave  and  Wind  Parameters 

Local  Variables: 

None. 

Subroutines  Called  from  LONG1  ( ):  None. 

LONG1  ( )  Called  from  Subroutines: 

MAIN  WAY 


Real 

Eddy  Viscosity  Coefficient 

Real 

Bottom  Friction  Coefficient 

Real 

Radiation  Stress  Coefficient 

Real 

Longshore  Wind  Stress  Coefficient 

Real 

Offshore  Water  Depth 

Real 

Roller  or  Bore  Term 

Integer 

Index  where  Wave  Probabilities 
Exceed  Threshold 

Integer 

Root  Mean  Square  Wave  Height 

Integer 

Wave  Number 
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Figure  36.  Subroutine  LONG1  Flowchart 
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5.36  Subroutine  MAIN  WAV 


Subroutine  Call: 

MAINWAV  (roller,  dxy,  xxl,  xshift,  b,  cl,  c2,  c3,  c4,  hrms,  xdelt_gr,  fqz,  nnn,  per,  xk,  fqd, 
Cg,  selfst,  theta,  theta2,  xtemp,  xktemp,  eblast,  htemp,  ptemp,  ebtemp,  iimax,  along,  blong,  clong, 
convg,  distmax,  rk,  bl,  surf,  j_ii,  dstart) 

Summary: 

Subroutine  MAIN_WAV  is  the  main  driver  for  coordinating  the  iterative  solution  method 
applied  to  solve  for  the  wave  and  current  parameters.  This  approach  is  necessary  because  several  of 
the  parameters  including  wave  height,  wave  length,  wave  celerity,  longshore  current  velocity,  and 
wave  induced  setup  are  interdependent,  as  well  as  depth  dependent. 


Input  Variables: 


b 

Real 

Empirical  Factor  in  Breaking  Model  =1.0 

cl 

Real 

Mixing/Eddy  Viscosity  Coefficient 

c2 

Real 

Bottom  Friction  Coefficient 

c3 

Real 

Factor  for  Radiation  Stress 

c4 

Real 

Friction  Coefficient  =  0.007 

Cg 

Real 

Wave  Group  Velocity 

dstart 

Real 

Starting  Depth  from  Input  File 

dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

fqd 

Real 

Peak  Frequency  at  the  Center  of  the 
Frequency  Band 

fqz 

Real 

Zero  Crossing  Frequency 

hrms 

Real 

Root  Mean  Square  Wave  Height 

nnn 

Integer 

Number  of  Points  in  Input  Depth  Array 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

self_st 

Char*  1 

Self  Start  Flag  (Yes  or  No) 

theta 

Real 

Wave  Angle 

xdelt  gr 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

xk 

Real 

Wave  Number 

xshift 

Real 

Horizontal  Cross-Shore  Location 

xxl  (points) 

Real 

Adjusted  Cross-Shore  Distances  from 
Depth  Profile 
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Output  Variables: 

along  (points) 

Real 

bl  (points) 

Real 

blong  (points) 

Real 

clong  (points) 

Real 

convg 

Logical 

distmax 

Real 

eblast 

Real 

ebtemp  (points) 

Real 

htemp  (points) 

Real 

iimax 

Integer 

jji 

Integer 

ptemp  (points) 

Real 

rk  (points, 4) 

Real 

surf 

Logical 

theta 

Real 

theta2 

Real 

xktemp  (points) 

Real 

xtemp  (points) 

Real 

Local  Variables: 

brklO 

Logical 

cg2 

Real 

check 

Real 

convcount 

Integer 

done 

Logical 

dp 

Real 

eb 

Real 

Horizontal  Mixing  Parameter 
Bottom  Slope 

Bottom  Friction  for  Deep  &  Shallow  Water 
Wind  Stress  Contribution  to 
Longshore  Current 
Energy  Equation  Convergence  Flag 
(True  or  False) 

Farthest  Offshore  Distance 
Roller  Dissipation  Term  at  Farthest 
Point  Offshore 

Temporary  Roller  Dissipation  Term 
Across  Transect 

Temporary  Variable  for  Significant  Wave 
Height  Values 

Number  of  Calculation  Locations 
Index  where  Wave  Probabilities 
Exceed  Threshold 
Percentage  of  Breaking  Waves  & 

Breaking  Types 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 
Flag  for  Low/No  Surf  Conditions 
(True  or  False) 

Wave  Angle 

Wave  Angle  at  Input  Starting  Depth 
Temporary  Variable  for  Wave  Number 
Temporary  Variable  for  Cross-Shore  Values 


Flag  Variable  to  Find  First  Location  Where  10% 

of  Waves  Are  Breaking  (True  or  False) 

Additional  Wave  Group  Velocity 

Difference  Between  Wave  Induced 

Setup  Calculations 

Number  of  Convergence  Iterations 

Loop  Control  Variable  for  Main  Wave 

Calculation  Loop  (True  or  False) 

Offshore  Water  Depth 
Temporary  Roller  Dissipation  Term 
Across  Transect 
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etanew  (points) 

Real 

etaold  (points) 

Real 

hrms2 

Real 

ii 

Integer 

11 

Real 

10 

Real 

P(4) 

Real 

pet (4) 

Real 

rhs 

Real 

thetaO 

Real 

xoff 

Real 

Subroutines  Called  from  MAIN_WAV  (  ): 

ABORT 

BALANCEQ 

GETBRK 

GETRHS 

LONG1 

PT_2 

SLFSTRT 

SETUP 

MAIN_WAV  ( )  Called  from  Subroutines: 

CALCSURF 


Wave  Induced  Setup  Estimated  at 
New  Location 

Wave  Induced  Setup  Estimated  at 
Previous  Location 

Wave  Height  for  Next  Onshore  Grid  Location 
Index  where  Wave  Probabilities 
Exceed  Threshold 

Wave  Length  at  Next  Onshore  Grid  Location 
Wave  Length  at  Grid  Cell  (1)  Offshore 
Array  for  Breaker  Percentage  Totals 
Percent  of  Different  Breaker  Types: 
pet  (1)  =  Spilling 
pet  (2)  =  Plunging 
pet  (3)  =  Surging 
pet  (4)  =  Total 

Right  Hand  Side  of  Energy  Balance  Equation 
Wave  Angle  at  Grid  (1)  Offshore 
Distance  Offshore 
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Figure  37.  Subroutine  MAIN_WAY  Flowchart 


i 
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5.37  Subroutine  MDSRF1 


Subroutine  Call: 


MDSRF1  (alfa,  chrlie,  pet,  echo,  foxtrt,  jgamma,  ihtll.  ihtl2.  file_out) 


Summary: 

Subroutine  MDSRF1  calculates  and  prints  the  modified  surf  index  number  to  the  output  file. 


Input  Variables: 

alfa 

Real 

Significant  Breaker  Height 

chrlie 

Real 

Dominant  Breaker  Period 

echo 

Real 

Breaker  Angle 

foxtrt 

Real 

Longshore  Current  Speed  and  Direction 

ihtll 

Real 

Wind  Speed 

ihtl2 

Real 

Wind  Direction 

jgamma 

Integer 

Temporary  Variable  set  to  Beach  Orientation 

pet (4) 

Real 

Percent  of  Different  Breaker  Types: 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

fileout 

Char*  40 

Output  File  Name 

Output  Variables: 

None. 

Local  Variables: 

idir 

Integer 

Index  for  Surf  Index  Wind  Direction 

index 

Integer 

Breaker  Type  Indicator  for  Surf  Index 

ispd 

Integer 

Index  for  Surf  Index  Wind  Speed  Lookup  in 
Modification  Table 

m 

Integer 

Temporary  Variable  to  Rotate  Direction 

percent 

Real 

Breaker  Type  Percentage 

srfmod 

Real 

Modified  Surf  Index  from  Sum  of  Values 
Resulting  from  Navy  Modification  Tables  in 
MDSRF2  ( ) 

sum 

Real 

Running  Total  of  Surf  Index 

suml 

Real 

Modified  Surf  Index  Value  for  Wave  Angle 

sum2 

Real 

Value  for  Longshore  Current 

temp 

Real 

Temporary  Wave  Angle  Variable 

theta2 

Real 

Rotated  Wind  Direction 

value 

Real 

Modification  Number 
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wind  (3,3,8) 


Real 


Surf  Index  Wind  Modification  Table 


Subroutines  Called  from  MDSRF1  ( ): 

MDSRF2 

MDSRF1  ( )  Called  from  Subroutines: 

SURF 
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Figure  38.  Subroutine  MDSRF1  Flowchart 
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5.38  Subroutine  MDSRF2 


Subroutine  Call: 

MDSRF2  (index,  xin.  yin,  value) 

Summary: 

Subroutine  MDSRF2  contains  the  modified  surf  index  (MSI)  tables.  The  MSI  number  is 
calculated  using  a  two  dimensional  linear  interpolation  by  areas. 

Input  Variables: 


index 

Integer 

Indicator  of  Breaker  Type 

xin 

Real 

X-Coordination  for  Surf  Index 
Modification  Matrix 

yin 

Real 

Y-Coordination  for  Surf  Index 
Modification  Matrix 

Output  Variables: 

value 

Real 

Returns  Modified  Surf  Index  Number 

Local  Variables: 

i 

Integer 

Loop  Counter  or  Array  Index 

il 

Integer 

Loop  Counter  or  .Array  Index 

i2 

Integer 

Loop  Counter  or  Array  Index 

ii 

Integer 

Loop  Counter  or  Array  Index 

ix  (4) 

Real 

All  Values  Set  to  1 1.00 

jy  (4) 

Real 

Values  Set  to  10.0,  11.0,  11.0,  9.0 

j 

Integer 

Loop  Counter  or  .Array  Index 

ji 

Integer 

Loop  Counter  or  Array  Index 

j2 

Integer 

Loop  Counter  or  Array  Index 

jj 

Integer 

Loop  Counter  or  .Array  Index 

tempi 

Real 

Temporary  Variable  Used  for  Interpolation 

x(ll) 

Real 

MSI  Indices 

xO  (4,11) 

Real 

Breaker  Period  Modification  table 

xdata 

Real 

Temporary  Index 

y  (ii) 

Real 

MSI  Indices 

yo  (4,i  i) 

Real 

Wave  Angle  Modification  table 

ydata 

Real 

Temporary  Index 

z  (11,11) 

Real 

Breaker  Modification  Matrix 

zO  (4,1 1,11) 

Real 

Whole  Breaker  Modification  Matrix 

zl  (40)-zl  1(40) 

Real 

Partial  Breaker  Modification  Arrays 
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zl2  (44) 
zzO  (484) 


Real 

Real 


Partial  Breaker  Modification  Array 
Equivalent  to  zO  (1,1,1)  zzO  (1) 


Subroutines  Called  from  MDSRF2  (  ):  None. 

MDSRF2  (  )  Called  from  Subroutines: 

MDSRF1 

Figure  39.  Subroutine  MDSRF2  Flowchart 


[  START  J 

I 

Setting  of  \ 
Surf  Index  \ 

Tables  / 

I 


Calculate 

X-Coordiante 
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5.39  Subroutine  NEW  BRK 
Subroutine  Call: 

NEW_BRK  (iimax,  bl,  rk,  htemp,  wid_ii,  p2) 

Summary: 

Subroutine  NEW_BRK  calculates  a  new  percentage  of  breaker  types  from  the  highest  10% 
of  the  wave  heights  (hrms)  when  the  bottom  slope  is  positive. 


Input  Variables: 


bl  (points) 

Real 

Bottom  Slope 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

iimax 

Integer 

Number  of  Calculation  Locations 

rk  (points,4) 

Real 

Matrix  of  Percentage  Breakers  and  Types 
Across  the  Transect 

widii 

Integer 

Offshore  Location  for  Surf  Zone  Width 

Output  Variables: 

p2  (4) 

Real 

Percent  of  Different  Breaker  Types  - 

Equivalent  to  pet  (4) 

p2  (1)  =  Spilling 

p2  (2)  =  Plunging 

p2  (3)  =  Surging 

p2  (4)  =  Total  ' 

Local  Variables: 

akl  (points) 

Real 

Temporary  Array  for  Wave  Height 

bkl  (points) 

Real 

Temporary  Array  Breaker  Type  =  1  Spilling 

bk2  (points) 

Real 

Temporary  Array  Breaker  Type  =  2  Plunging 

bk3  (points) 

Real 

Temporary  Array  Breaker  Type  =  3  Surging 

bk4  (points) 

Real 

Temporary  Array  Breaker  Type  =  4  Total 

Total  Percentage  of  Breakers 

i 

Integer 

Loop  Counter 

ii 

Integer 

Loop  Counter 

nval 

Integer 

Number  of  Positive  Slope  Occurrences 

xl 

Real 

0.1  %  of  Highest  Breakers  to  Examine  for  Type 
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x2 


Integer 


Loop  Limit  -  Set  to  Top  Percentage  of  Significant 
Wave  Height  Values 


Subroutines  Called  from  NEWBRK  ( ): 

ABORT 

RN2 

NEW  BRK  (  )  Called  from  Subroutines: 

SHORTOUT 
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Figure  40.  Subroutine  NEW  BRK  Flowchart 


START 

T 


Number  of 
Coordinate 
Positions  Loop 


I 


Number  of 
HRMS  Valves  in 
Top  10  %  Loop 


l 

Calculate  Percent  Spilling 


\ 

Calculate  Percent  Plunging 


1 

Calculate  Percent  Surging 


I 


Calculate  Total  Percentage 
of  Breakers 


i 


End  of  Loop 


I 

Set  Number  of 
HRMS  Values  Within 
the  Top  1 0  % 

I 

Initialize  Wave  Type 
Percentage  Array 


\ 


Set  all  Breaker  Type 
Summation  Percentages 


I 

(  RETURN  ] 
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5.40  Subroutine  NONLIN 


Subroutine  Call: 

NONLIN  (j_ii,  xktemp,  htemp,  dxy,  ebtemp,  thetaZ  cf,  iimax,  v) 

Summary: 

Subroutine  NONLIN  calculates  the  cross-shore  distribution  of  the  longshore  current  using  a 
nonlinear  bed  stress  as  the  restoring  force  in  the  momentum  equation. 


Input  Variables: 


cf 

Real 

Coefficient  of  Friction  for  the  Bottom  Stress 

dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

ebtemp  (points) 

Real 

Roller  Dissipation  Term  Across  Transect 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

iimax 

Integer 

Number  of  Calculation  Locations 

jji 

Integer 

Index  where  Wave  Probabilities 

Exceed  Threshold 

theta2 

Real 

Wave  Angle  at  Input  Starting  Depth 

xktemp  (points) 

Real 

Array  for  Wave  Number 

Output  Variables: 

v  (points) 

Real 

Longshore  Current  Velocity  Distribution 

Local  Variables: 

c 

Real 

Temporary  Variable  Used  for  Longshore  Current 
Calculation 

c3 

Real 

Refraction  Coefficient  Based  upon  Farthest 
Offshore  Wave  Angle 

dp 

Real 

Offshore  Water  Depth 

ebn 

Real 

Temporary  Roller  Dissipation  Term 

grdjpt 

Integer 

Loop  Counter 

hrms 

Real 

Root  Mean  Square  Wave  Height 

q 

Real 

Longshore  Current  Momentum  Flux 

vtmp 

Real 

Temporary  Longshore  Current  Velocity 

xk 

Real 

Wave  Number 
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Subroutines  Called  from  NONLIN  (  ): 

NONLIN2 

NONLIN  ( )  Called  from  Subroutines: 

CALC  SURF 

Figure  41.  Subroutine  NONLIN  Flowchart 


Calculate  Longshore 
Current  Values 


Call 

NONLIN2 


\ 


Calculate  Longshore 
Current  Velocity 

T 


End  of  Loop 


(  RETURN  ] 


5.41  Subroutine  NONLIN2 


Subroutine  Call: 

N0NLIN2  (xk,  hrras,  dp,  theta2,  q,  v) 

Summary: 

Subroutine  NONLEN'2  initializes  variables  in  the  longshore  momentum  equation  and  checks 
for  convergence  of  the  iterative  solution  method. 


Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

turns 

Real 

Root  Mean  Square  Wave  Height 

q 

Real 

Longshore  Current  Momentum  Flux 

theta2 

Real 

Wave  Angle  at  Input  Starting  Depth 

xk 

Real 

Wave  Number 

Output  Variables: 

v  Real  Longshore  Current  Velocity 

Local  Variables: 


convg 

Logical 

Convergence  Flag  (True  or  False) 

freq 

Real 

Wave  Frequency 

h 

Real 

Wave  Height 

kount 

Integer 

Counter 

1 

Real 

Wave  Length 

t 

Real 

Wave  Period 

u 

Real 

Mean  Cross-Shore  Current  Velocity 

Subroutines  Called  from  NONLIN2  (  ): 

ABORT 

CONANG 

NONLIN2  ( )  Called  from  Subroutines: 
NONLIN 
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Figure  42.  Subroutine  N0NLIN2  Flowchart 
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5.42  Subroutine  PERCENT 


Subroutine  Call: 

PERCENT  (hmis,  period,  dp,  slope,  p) 

Summary: 

Subroutine  PERCENT  calculates  the  percentage  of  each  type  of  breaking  wave  in  the  surf 


zone. 

Input  Variables: 


dp 

Real 

hrms 

Real 

period 

Real 

slope 

Real 

Output  Variables: 

P(4) 

Real 

Local  Variables: 

ffac  (3) 

Real 

gtemp 

Real 

hhigh 

Real 

hlow 

Real 

integrat 

Real 

pflag 

Logical 

param 

Real 

Offshore  Water  Depth 
Root  Mean  Square  Wave  Height 
Peak  Period 
Bottom  Slope 


Array  of  Percentage  of  Breaker  Types 

pet  (1)  -  Spilling 

pet  (2)  -  Plunging 

pet  (3)  -  Surging 

pet  (4)  -  Total  Percentage 


Array  for  Percentage  Breaker  Totals 
Gravity 

Upper  Bound  of  Integration 
Lower  Bound  of  Integration 
Wave  Height  Distribution  Calculated  at  a 
Specific  Location 

Weighting  Factor  Flag  (True  or  False) 
Integral  Multiplier 


Subroutines  Called  from  PERCENT  (  ): 

GET  P 
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Functions  Called  from  PERCENT  (  ): 

INTEGRA! 


PERCENT  ( )  Called  from  Subroutines: 

GETBRK 

SLFSTRT 

Figure  43.  Subroutine  PERCENT  Flowchart 


[  START  J 

I 

<Scmng 
Integral 
Multiplier 


I 


Page  106 


5.43  Subroutine  PRT  OUT1 


Subroutine  Call: 

PRT_OUTl  (j_ii,  xdelt.  iimax,  dxy,  xtemp,  xktemp,  htemp,  ptemp,  v) 

Summary: 

Subroutine  PRT  OUT1  prints  columnar  data,  cross-shore  distributions  of  wave  and  surf 
parameters,  to  the  detailed  SURF  3.1  output  file  when  requested  by  the  user.  This  data  is  interpolated 
to  the  user  defined  grid  step,  if  possible. 


Input  Variables: 


dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

j_ii 

Integer 

Index  where  Wave  Probabilities 

Exceed  Threshold 

iimax 

Integer 

Number  of  Calculation  Locations 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

ptemp  (points) 

Real 

Percentage  of  Breaking  Waves  & 

Breaker  Types 

v  (points) 

Real 

Longshore  Current  Velocity 

xdelt 

Real 

Surf  Zone  Output  Interval 

xktemp  (points) 

Real 

Temporary  Wave  Number  Array 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

Output  Variables: 

None. 

Local  Variables: 

dpi 

Real 

Offshore  Depth 

hmax 

Real 

Maximum  Wave  Height 

houtl 

Real 

Significant  Wave  Height 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 

ii 

Integer 

Array  Index  Number 

jj 

Integer 

Iteration  Count 

pbreak 

Real 

Percentage  Breaking  Waves 

vlngl 

Real 

Longshore  Current  Velocity 

wlen 

Real 

Wave  Length 

xoffl 

Real 

Distance  Offshore 
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Subroutines  Called  from  PRT  OUT1  ( ): 

ABORT 

GT_P 

GRIDOUT 

PRT_OUTl  ( )  Called  from  Subroutines: 

CALCSURF 

Figure  44.  Subroutine  PRT  OUT1  Flowchart 
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5.44  Subroutine  PRT  OUT2 
Subroutine  Call: 

PRT  OUT2  (j_ii,  xdelt,  iimax,  dxy,  xtemp,  xktemp,  htemp,  ptemp,  v) 

Summary: 

Subroutine  PRTOUT2  writes  the  detailed  surf  output. 


Input  Variables: 

dxy  (points) 

Real 

jji 

Integer 

iimax 

Integer 

htemp  (points) 

Real 

ptemp  (points) 

Real 

v  (points) 

Real 

xdelt 

Real 

xktemp  (points) 

Real 

xtemp  (points) 

Real 

Output  Variables: 

None. 

Local  Variables: 

dpi 

Real 

hmax 

Real 

houtl 

Real 

hrmsl 

Real 

ii 

Integer 

jj 

Integer 

pbreak 

Real 

vlng 

Real 

wlen 

Real 

xoffl 

Real 

Corresponding  Depths  with  Tide 
Index  where  Wave  Probabilities 
Exceed  Threshold 
Number  of  Calculation  Locations 
Temporary  Variable  for  Significant  Wave 
Height  Values 

Percentage  of  Breaking  Waves  and 

Breaker  Types 

Longshore  Current  Velocity 

Surf  Zone  Output  Interval 

Temporary  Wave  Number  Array 

Temporary  Variable  for  Cross-Shore  Values 


Offshore  Depth 
Maximum  Wave  Height 
Significant  Wave  Height 
Root  Mean  Square  Wave  Height 
Array  Index  Number 
Iteration  Counter 
Percentage  Breaking  Waves 
Longshore  Current  Velocity 
Wave  Length 
Distance  Offshore 


Subroutines  Called  from  PRT  OUT2  ( ): 

GTJP 

ABORT 
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PRT  OUT2  ( )  Called  from  Subroutines: 

CALCSURF 

Figure  45.  Subroutine  PRT  OUT2  Flowchart 
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5.45  Subroutine  PRT  OUT3 
Subroutine  Call: 

PRT  OUT3  (file_dat) 

Summary: 

Subroutine  PRT  OUT3  writes  out  the  detailed  output  from  the  model. 

Input  Variables: 

file  dat  Char*40  Output  File  name  *.dat 

Output  Variables:  None. 

Local  Variables: 

line  Char*80  Temporary  String 

Subroutines  Called  from  PRT  OUT3  ( ):  None. 

PRT  OUT3  ( )  Called  from  Subroutines: 

SURF 

Figure  46.  Subroutine  PRT  OUT3  Flowchart 
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5.46  Subroutine  PT2 


Subroutine  Call: 

PT2  (10.  thetaO,  fqd,  dp,  theta,  xk,  1,  Cg) 

Summary: 

Subroutine  PT2  calculates  wave  parameters  from  linear  theory  relations. 


Cg  -  nC 


Group  Velocity 


n  =  {  [1  + 


2_kh_ 
sinh  2  k  h 


sin  9  _  sin  Q0 

C  Co  Wave  angle  from  Snell’s  law 


Input  Variables: 

dp 

Real 

Offshore  Water  Depth 

fqd 

Real 

Peak  Frequency 

10 

Real 

Wave  Length  at  Offshore  Point 

thetaO 

Real 

Wave  Angle  at  Offshore  Point 

xk 

Real 

Wave  Number 

Output  Variables: 

Cg 

Real 

Group  Velocity 

1 

Real 

Wave  Length 

theta 

Real 

Wave  Angle 

xk 

Real 

Wave  Number 

Local  Variables: 

c 


Real 


Temporary  Variable 


Subroutines  Called  from  PT2  ( ): 

WAVENUM 


PT2  (  )  Called  from  Subroutines: 
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MAIN_WAV 

SLFSTRT 

Figure  47.  Subroutine  PT2  Flowchart 
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5.47  Subroutine  RAD  ST1 


Subroutine  Call: 

RAD_ST1  (ifreq,  freq,  idirec,  xffom,  esowm,  freql,  ffeq2,  dstart,  igamma,  theta,  hrms,  surf, 
fqd,  per,  fqz) 

Summary: 

Subroutine  RAD_ST1  searches  the  directional  wave  spectrum  to  identify  the  dominant  wave 
frequency  and  sums  the  wave  energy  directed  toward  shore.  The  flux  of  momentum  or  Radiation 
Stress,  which  contributes  to  driving  the  longshore  current,  is  calculated  following  Thornton  and  Guza 


SXy  (0,j)  =  E  (6, fin  (f)  sin  a  (f)  cos  a  (f) 


(1986). 

In  the  above  equation  Sxy  is  the  Radiation  Stress,  E  is  the  total  energy  in  the  directional  wave 
spectrum,  n  is  the  ratio  of  wave  group  velocity  to  wave  velocity,  and  a  is  the  wave  angle.  The  ratio 


n 


Q 

C 


0.5  (1  + 


2  k  h 
sinh  kh 


n  from  linear  wave  theory  is  given  by: 

where,  Cg  is  the  group  velocity,  C  is  the  wave  velocity  or  celerity,  k  is  the  wave  number  and  h  is  the 
local  water  depth. 

Input  Variables: 

dstart  Real 

esowm  (dirNum,ffeqNum)  Real 
ffeq  (ffeqNum)  Real 

freql(freqNum)  Real 


Input  Starting  Depth 

Directional  Wave  Spectrum 

Input  Wave  Spectrum  Center  Frequencies 

Beginning  Frequency  Bin  Values 
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freq2  (ffeqNum) 

Real 

idirec 

Integer 

ifreq 

Integer 

igamma 

Integer 

xffom  (freqNum) 

Real 

Output  Variables: 

fqd 

Real 

fqz 

Real 

hrms 

Real 

per 

Real 

surf 

Logical 

theta 

Real 

Local  Variables: 

direc 

Real 

ees 

Real 

esum 

Real 

esumm 

Real 

frd 

Real 

idir 

Integer 

ifrq 

Integer 

m 

Integer 

maxfrq 

Integer 

summax 

Real 

sumzero 

Real 

sxy 

Real 

sxysum 

Real 

temp 

Real 

temp2 

Real 

theta2 

Real 

xk 

Real 

xkd 

Real 

Ending  Frequency  Bin  Values 
Number  of  Directions  in  Input  Spectrum 
Number  of  Frequencies  in  Input  Spectrum 
Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 
Direction  Array,  Direction  Wave  Energy  Comes 
From 


Peak  Frequency  at  the  Center  of  the 

Frequency  Band 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Peak  Period  of  Directional  Wave  Spectrum 

Flag  for  Low  or  No  Surf  Conditions 

(True  or  False) 

Wave  Angle 


Wave  Direction 

Spectral  Density  at  a  Particular 

Frequency  and  Direction 

Sum  of  Energy  in  One  Frequency  Band 

Over  all  Directions 

Sum  of  All  Energy  in  Directional  Spectrum 
Wave  Frequency 
Loop  Counter 
Loop  Counter 

Temporary  Variable  for  Rotating  Wave  Angle 
Frequency  at  Maximum  Spectral  Density 
Frequency  Band  with  Maximum  Energy 
Summation  of  Zero-Crossing  Frequency  Energy 
Radiation  Stress 
Sum  of  Radiation  Stress  Energy 
Temporary  Variable  in  Radiation 
Stress  Calculation 

Temporary  Variable  for  Frequency  Band  with 
Maximum  Energy 

Angle  Between  Wave  Ray  and  Beach 
Perpendicular  Projection 
Wave  Number 

Wave  Number  Multiplied  by  the 
Local  Water  Depth 
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Subroutines  Called  from  RAD  ST1  (  ): 

RADST2 

WAVENUM 

RAD_ST1  ( )  Called  from  Subroutines: 

CALCSURF 
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Figure  48.  Subroutine  RAD_ST1  Flowchart 
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5.48  Subroutine  RAD  ST2 
Subroutine  Call: 

RAD_ST2  (freq,  sxysum,  sumzero,  esumm,  maxfrq,  dstart,  theta,  hrms,  surf,  fqd,  per,  fqz) 

Summary: 

Subroutine  RAD_ST2  calculates  several  parameters  based  on  the  total  energy  in  the 
directional  wave  spectrum.  A  check  is  performed  to  confirm  that  wave  energy  is  directed  onshore 
before  writing  summary  information  to  the  output  file. 


Input  Variables: 

dstart 

Real 

esumm 

Real 

freq  (freqNum) 

Real 

maxfrq 

Integer 

sumzero 

Real 

sxysum 

Real 

Output  Variables: 

fqd 

Real 

fqz 

Real 

hrms 

Real 

per 

Real 

surf 

Logical 

theta 

Real 

Local  Variables: 

hs 

Real 

sxy2 

Real 

temp 

Real 

theta3 

Real 

xk 

Real 

xkd 

Real 

Input  Starting  Depth 

Sum  of  All  Energy  in  Directional  Spectrum 
Input  Wave  Spectrum  Center  Frequencies 
Frequency  at  Maximum  Spectral  Density 
Summation  of  Zero-Crossing 
Frequency  Energy 
Sum  of  Radiation  Stress  energy 


Peak  Frequency 

Zero  Crossing  Frequency 

Root  Mean  Square  Wave  Height 

Peak  Period  of  Directional  Wave  Spectrum 

Logical  Flag  for  Low/No  Surf  Conditions 

(True  or  False) 

Wave  Angle 


Significant  Wave  Height 

Temporary  Wave  Energy 

Temporary  Variable  for  Energy 

Wave  Angle  in  Degrees 

Wave  Number  Calculated  at  Peak  Frequency  and 

Input  Starting  Depth 

Wave  Number  *  Water  Depth 


Page  1 1 8 


Subroutines  Called  from  RAD_ST2  ( ): 

S_NOSURF 

WAVENUM 

RAD  ST2  ( )  Called  from  Subroutines: 

RADST1 

Figure  49.  Subroutine  RAD  ST2  Flowchart 
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5.49  Subroutine  READRFRC 


Subroutine  Call: 

READRFRC  (fracname,  idwsfreq,  idwsdirec,  xcoefF,  xtheta,  sfreq,  sdir) 

Summary: 

Subroutine  READRFRC  reads  refraction  information  from  a  formatted  input  file.  The 
matrices  contained  in  these  files  are  used  to  shoal  and  refract  a  directional  wave  spectrum  from  an 
offshore  point  to  a  location  where  depth  information  is  available.  The  number  of  frequency  bins  must 
not  exceed  50  and  the  number  of  direction  bins  must  not  exceed  ISO.  The  directional  coverage  of  the 
refraction  and  shoaling  coefficients  must  range  from  0  to  360  degrees.  Partial  coverage  over  a 
fraction  of  the  compass  (e.g.  180  degree  sector)  will  introduce  errors. 

Input  Variables: 

fracname  Char*40  Wave  Refraction  File 


Output  Variables: 


idwsdirec 

Integer 

Number  of  rows  (Directions)  in  the 
Directional  Wave  Spectrum  Matrix 

idwsfreq 

Integer 

Number  of  columns  (Frequencies)  in  the 
Directional  Wave  Spectrum  Matrix 

sdir  (dirNum) 

Real 

Direction  Array  for  each  bin  in  the 
Directional  Wave  Spectrum 

sfreq  (freqNum) 

Real 

Center  Frequency  of  each  Directional 
Wave  Spectrum 

xcoeff  (dirNum, freqNum) 

Real 

Wave  Height  Refraction  Coefficients 

xtheta  (dirNum, freqNum) 

Real 

Angle  Refraction  Coefficients 

Local  Variables: 

cfmatch 

Logical 

Flag  for  Center  Frequency  Match 

cffeq  (freqNum) 

Real 

Center  Frequency  of  each  Bin 

col 

Real 

Number  of  Columns 

dangle 

Real 

Angle  Between  Directional  Bins 

dir 

Real 

Number  of  Angles 

dirin 

Integer 

X-Coordinates  of  known  values 

dirord 

Integer 

Direction  of  Waves 
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dirouts  (dirNum) 

Real 

dirs  (dirNum) 

Real 

dmatch 

Logical 

dots 

Integer 

drl 

Real 

dth 

Real 

dum 

Real 

dumstr 

Char*  80 

fmatch 

Logical 

fnum 

Integer 

found 

Integer 

frchk 

Integer 

frq 

Real 

I 

Integer 

ii 

Integer 

icol 

Integer 

idir 

Integer 

idirec 

Integer 

ifreq 

Integer 

ifrq 

Integer 

instat 

Integer 

irow 

Integer 

j 

Integer 

jj 

Integer 

k 

Integer 

kk 

Integer 

Ifreq 

Real 

lowcut 

Integer 

mpnt 

Integer 

refs  (dirNum) 

Real 

rfrtmp  (dirNum,ffeqNum) 

Real 

row 

Real 

rtmpout  (dirNum) 

Real 

sfreqin  (dirNum) 

Real 

shltmp  (dirNum,freqNum) 

Real 

splout  (dirNum) 

Real 

stmpout  (dirNum) 

Real 

temp  (dirNum, ffeqNum) 

Real 

temp2  (dirNum,ffeqNum) 

Real 

tmpinr  (dirNum) 

Real 

1  -  Direction  Waves  are  coming  from 

2  -  Direction  Waves  are  going  to 
Interpolated  X-Coordinates 

Temporary  Direction  Wave  Energy  Comes  From 
Flag  for  Directional  Match 
Y-Coordinates  of  known  values 
Initial  Direction  Bin 

Temporary'  Angle  Between  Directional  Bins 

Tempo  rary  Variable 

Temporary  Variable 

Flag  for  Frequency  Match 

Bin  Number 

Flag  Indicator 

Total  Number  of  Frequencies 
Number  of  Frequencies 
Loop  Counter 
Counter 

Number  of  Columns 
Loop  Counter 

Number  of  Rows  (Directions)  in  the 

Refraction  Shoaling  Matrix 

Number  of  Columns  (Frequencies)  in  the 

Refraction  Shoaling  Matrix 

Loop  Counter 

Error  S  tarns 

Number  of  Rows 

Loop  Counter 

Counter 

Counter 

Counter 

Lower  Frequency  Bin  Limit 
Lower  Cut  Off  Limit 
Number  of  Rows  divided  by  2 
Temporary'  Array 

Temporary'  Matrix  for  Reversing  Wave 

Direction 

Number  of  Rows 

Interpolated  Coordinates 

Temporary  Frequency  Array 

Temporary  Matrix  for  Reversing  Wave 

Direction 

interpolated  Y-Coordinates 
Interpolated  Coordinates 
Temporary  Variable 
Temporary  Variable 
Temporary  Variable 
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tmpins  (dirNum) 

Real 

Temporary  Variable 

ufreq 

Real 

Upper  Frequency  Bin  Limit 

upcut 

Integer 

Upper  Cut  Off  Limit 

xfrom  (dirNum) 

Real 

Direction  Wave  Energy  Comes  From 

Subroutines  Called  from  READRFRC  ( ): 


ABORT 

GENRLSPL 

READRFRC  ( )  Called  from  Subroutines: 

SURF 
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Figure  50.  Subroutine  READRFRC  Flowchart 
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5.50  Subroutine  READSPEC 


Subroutine  Call: 

READSPEC  (ifreq,  idirec,  Cfreq.  Lfreq,  Ufreq,  xfrom,esowm,  period,  ehsig,  dangle,  spefile) 

Summary: 

Subroutine  READSPEC  opens  and  reads  a  directional  wave  spectrum  file,  which  must 
conform  to  a  specific  format,  but  the  number  of  frequencies  and  directions  can  vary.  The  maximum 
number  of  directions  is  180  and  the  maximum  number  of  frequencies  is  50.  The  directions  should 
be  evenly  spaced,  and  the  frequency  bins  can  be  fixed  or  variable  width  with  units  of  energy  density 
(mA2/(Hz*radians)).  This  energy  density  matrix  is  initialized,  filled,  and  converted  to  units  of  feet 
squared  inside  this  subroutine.  Also,  the  direction  of  wave  energy  can  be  the  direction  FROM  which 
waves  are  coming  or  TO  which  waves  are  going  as  denoted  in  the  tenth  header  line  by  a  1  or  2 
respectively.  The  directional  wave  spectrum  must  be  defined  from  0  to  360  degrees.  Use  of  partial 
directional  sectors  (e.g.  0  to  180  degrees)  will  cause  errors. 


Input  Variables: 

None. 

Output  Variables: 

Cfreq  (freqNum) 

Real 

dangle 

Real 

ehsig 

Real 

esowm  (dirNum, freqNum) 

Real 

idirec 

Integer 

ifreq 

Integer 

Lfreq  (freqNum) 

Real 

period  (freqNum) 

Real 

spefile 

Char*40 

Ufreq  (freqNum) 

Real 

xffom  (dirNum) 

Real 

Center  Frequency  Bin  Limit 

Angle  Between  Directional  Bins 

Significant  Wave  Height  from 

Directional  Spectrum 

Directional  Wave  Spectrum 

Number  of  Direction  Bins  in  Input  Spectrum 

Number  of  Frequency  Bins  in  Input  Spectrum 

Lower  Frequency  Bin  Limit 

Period  Array  (1 /Frequency) 

Wave  Spectrum  File  Name 
Upper  Frequency  Bin  Limit 
Direction  Array,  Direction  Wave  Energy 
Comes  From 
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Local  Variables: 


col 

Real 

df 

Real 

dir 

Real 

dirord 

Integer 

dth 

Real 

dum 

Char*l 

drl 

Real 

fnum 

Integer 

frq 

Real 

ftsq2msg 

Real 

I 

Integer 

icol 

Integer 

idir 

Integer 

ifrq 

Integer 

instat 

Integer 

irow 

Integer 

j 

Integer 

mpnt 

Integer 

mult 

Real 

row 

Real 

temp  (dirNum.dirNum) 

Real 

Number  of  Columns 

Difference  between  Upper  &  Lower  Bins 
Number  of  Angles 
Direction  of  Waves 

1  -  Direction  Waves  are  coming  from 

2  -  Direction  Waves  are  going  to 
Width  of  Direction  Bin 
Temporary  Variable 

Initial  Direction  Bin 
Bin  Number 
Number  of  Frequencies 
Conversion  Factor 
Loop  Counter 
Number  of  Columns 
Direction  Loop  Counter 
Loop  Counter 
Error  Status 
Number  of  Rows 
Loop  Counter 

Number  of  Rows  divided  by  2 
Temporary  Calculation  Variable 
Number  of  Rows 
Temporary  Array 


Subroutines  Called  from  READSPEC  ( ): 

ABORT 

GTSIGH 

READSPEC  (  )  Called  from  Subroutines: 

SURF 


Page  126 


Figure  51.  Subroutine  READSPEC  Flowchart 


Flip  Directoinal 
Wave  Spectrum 
Matrix  180° 
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5.51  Subroutine  REFRAC 
Subroutine  Call: 

REFRAC  (idirec,  ifreq.  dangle,  xtheta,  xcoeff,  esowrn,  ehsig) 

Summary: 

For  each  frequency  and  direction  bin  in  the  input  directional  wave  spectrum,  the  shallow  water 
direction  band  for  each  deep  water  direction  band  is  found.  Wave  energy  from  each  deep  water  band 
is  multiplied  by  the  combined  refraction/shoaling  coefficient  and  moved  into  the  proper  shallow  water 
band  to  provide  a  shallow  water  directional  spectrum. 

Input  Variables: 

dangle  Real 

idirec  Integer 

ifreq  Real 

xcoeff  (dirNum,freqNum)  Real 

xtheta  (dirNum,freqNum)  Real 

Output  Variables: 

ehsig  Real  Significant  Wave  Height  from 

Directional  Spectrum 

esowrn  (dirNum,ffeqNum)  Real  Directional  Wave  Spectrum 

Local  Variables: 


Angle  Between  Directional  Bins 
Number  of  Direction  Bins  in  Input  Spectrum 
Number  of  Frequencies  in  Input  Spectrum 
Wave  Height  Refraction  Coefficients 
Angle  Refraction  Coefficients 


esite  (dirNum,freqNum)  Real 

idir  Integer 

ifrq  Integer 

itemp  Integer 

itheta  (dirNum,freqNum)  Integer 
jdir  Integer 

mtemp  Integer 

sum  Real 

sum2  Real 

ytheta  Real 


Subroutines  Called  from  REFRAC  ( ): 


Directional  Spectrum  in  Shallow  Water 
Direction  Loop  Counter 
Frequency  Loop  Counter 
Temporary  Wave  Angle  Variable 
Shoreward  Energy  Spectrum 
Loop  Variable 

Temporaiy  Wave  Angle  Variable 
Temporary  Wave  Energy  Summation  Variable 
Temporary  Wave  Energy  Summation  Variable 
Temporary  Wave  Angle  Variable 

None. 
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REFRAC  ( )  Called  from  Subroutines: 

SURF 

Figure  52.  Subroutine  REFRAC  Flowchart 
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5.52  Subroutine  RN2 


Subroutine  Call: 

RN2  (n,  x,  yl,  y2,  y3,  y4) 

Summary: 

Subroutine  RN2  calculates  percentages  of  each  type  of  breaker  in  the  surf  zone. 

Input  Variables: 


n 

Integer 

Number  of  Waves  Considered  Breaking  on  £ 
Positive  Bottom  Slope 

x  (points) 

Real 

Temporary  Significant  Wave  Height  Array 

yl  (points) 

Real 

Spilling  Breaker  Type 

y2  (points) 

Real 

Plunging  Breaker  Type 

y3  (points) 

Real 

Surging  Breaker  Type 

y4  (points) 

Real 

Total  Number  of  Breakers 

Output  Variables: 

yl  (points) 

Real 

Spilling  Array  Breaker  Type 

yl  (points) 

Real 

Plunging  Array  Breaker  Type 

y3  (points) 

Real 

Surging  Array  Breaker  Type 

y4  (points) 

Real 

Total  Array  Breaker  Type 

Local  Variables: 

hold 

Real 

Temporary  Variable  Used  for  Repositioning 

i 

Integer 

Loop  Counter 

j 

Integer 

Loop  Counter 

js 

Integer 

Loop  Starting  Index 

m 

Integer 

Number  of  Waves  Considered  Breaking  on  a 
Positive  Slope 

Subroutines  Called  from  RN2  ( ): 

None. 

RN2  (  )  Called  from  Subroutines: 

NEWJBRK 
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5.53  Subroutine  S  COEFF 


Subroutine  Call: 

S_COEFF  (dp,  fqd,  hrms,  theta,  c,  xk,  wdir,  igamma,  wdspd,  cl,c2,  c3,  c4,  cf,  vwind) 

Summary: 

Subroutine  S_COEFF  calculates  several  parameters  in  the  longshore  current  equation 
including  the  Radiation  Stress,  the  bottom  stress,  and  the  wind  stress.  A  check  is  performed  to  assure 
that  wave  induced  motion  is  not  dominated  by  wind  effects  and  a  warning  message  is  written  to  the 
output  file  if  this  condition  is  violated.  An  assumption  is  made  that  if  the  wave  induced  orbital 
velocity  is  greater  than  the  wind-forced  component  of  the  longshore  current,  the  local  conditions  are 
wave  dominated. 


Input  Variables: 


c 

Real 

Wave  Celerity  at  Input  Starting  Depth 

dp 

Real 

Water  Depth  Offshore 

fqd 

Real 

Peak  Frequency  from  Directional  Spectrum 

hrms 

Real 

Root  Mean  Square  Wave  Height 

igamma 

Integer 

Beach  Orientation,  Compass  Heading 
Directly  Toward  Beach 

theta 

Real 

Wave  Angle 

wdir 

Real 

Input  Wind  Direction  Compass  Heading 

wdspd 

Real 

Input  Wind  Speed 

xk 

Real 

Wave  Length  at  Input  Starting  Depth 

Output  Variables: 

cl 

Real 

Mixing/Eddy  Viscosity  Coefficient 

c2 

Real 

Bottom  Friction  Coefficient 

c3 

Real 

Factor  for  Radiation  Stress 

c4 

Real 

Friction  Coefficient 

vwind 

Real 

Wind  Driven  Longshore  Current  Velocity 
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Local  Variables: 


c4tmp 

Real 

Temporary  Variable  Used  in  Wind  Velocity 
Vector  Calculation 

cd 

Real 

Coefficient  of  Drag  Used  in  Wind 

Velocity  Calculation 

cf 

Real 

Coefficient  of  Friction  for  the  Bottom 

dwind 

Real 

Sign  of  Wind  Vector  (Positive  or  Negative) 

m 

Integer 

Temporary  Variable  Used  in  Rotating 

Wind  Angle 

theta4 

Real 

Rotated  Wind  Direction 

uorb 

Real 

Wave  Particle  Orbital  Velocity 

xn 

Real 

Eddy  Viscosity  Mixing  Coefficient 

Subroutines  Called  from  S 

COEFF  ( ): 

None. 

S  COEFF  (  )  Called  from  Subroutines: 

CALCSURF 
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Figure  54.  Subroutine  S  COEFF  Flowchart 


I 


[Calculate  Drag  Coefficient] 
for  Wind  Velocity 


\ 
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I 
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^article  Orbital  Velocity? 


I 


Yes 


Warning  Message: 
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I 
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5.54  Subroutine  SNOSURF 
Subroutine  Call: 

S_NOSURF  ( hsig,  surf) 

Summary: 

Subroutine  S_NOSURF  is  called  to  determine  if  local  conditions  are  significant  enough  to 
proceed  with  surf  zone  calculations.  The  minimum  condition  for  continuation  is  that  the  significant 
wave  height  calculated  from  the  directional  wave  spectrum  must  be  greater  than  0.15  m. 


Input  Variables: 

hsig 

Real 

Significant  Wave  Height 

Output  Variables: 

surf 

Logical 

Flag  to  Indicate  Low  or  No  Surf  Conditions 
(True  or  False) 

Local  Variables: 

None. 

Subroutines  Called  from  S_NOSURF  ( ):  None. 

S  NOSURF  (  )  Called  from  Subroutines: 

CALCSURF 
RAD  ST2 
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5.55  Subroutine  S  TIDE 


Subroutine  Call: 

STIDE  (tide,  ydepth,  nnn,  dxyl,  xxl,  dxy,  xshift) 

Summary: 

Subroutine  S  TIDE  adds  the  tidal  elevation  to  each  cross-shore  point  in  the  input  depth 

profile. 

Input  Variables: 

dxyl  (points) 
nnn 
tide 

xxl  (points) 
ydepth 

Output  Variables: 

dxy  (points)  Real  Adjusted  Depths  with  Tide 

xshift  Real  Offshore  Distance 

Local  Variables: 

ddiff  Real  Change  in  Water  Depth 

n  Integer  Loop  Counter 

nn  Integer  Loop  Counter 

mm  Integer  Loop  Counter 

xdiff  Real  Change  in  Cross-Shore  Location 

ztide  Real  Tide  Level 

Subroutines  Called  from  S  TIDE  ( ):  None. 

S  TIDE  ( )  Called  from  Subroutines: 

CALCSURF 


Real 

Corresponding  Depths  without  Tide 

Integer 

Number  of  Points  in  Input  Depth  Array 

Real 

Tide  Level 

Real 

Adjusted  Cross-Shore  Distances  from 

Depth  Profile 

Char*  1 

Usage  of  Input  Depth  (Yes/No) 
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Figure  56.  Subroutine  S  TIDE  Flowchart 
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5.56  Subroutine  SEAFIT 


Subroutine  Call: 

SEAFIT  (hsig,  per,  dir,  ifreq,  idirec,  freql,  freq2,  xfrom,  esowm) 

Summary: 

Subroutine  SEAFIT  calculates  a  directional  wave  spectrum  from  an  input  wave  height  and 
wave  period  using  a  Pierson-Moskowitz  spectrum  representation  and  a  cosine  to  the  fourth  directional 
spreading  function.  The  modified  Pierson-Moskowitz  equation  (from  Pierson  and  Moskowitz,  1964) 


E(f)  =  ag2  w5 


g[-bt  w0/ H-/y 


provides  wave  energy  at  each  frequency  from  the  following  equation: 
where  : 


W  =  27lf 


a  =  0.0081 

b  =  0.74 

_  £ 
w0-- 

in  which/is  the  wave  frequency  in  Hertz,  g  is  gravity,  and  U  is  the  wind  speed  in  meters  per  second 
measured  at  19.5  m  above  the  sea  surface.  The  spectrum  E(f)  is  a  vector  of  spectral  densities  and  it 
is  assumed  that  each  density  is  integrated  from  the  lower  limit  of  the  frequency  bin  to  the  upper  limit 
of  the  frequency  bin. 
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Input  Variables: 


dir 

Real 

Wave  Direction 

freql  (freqNum) 

Real 

Beginning  Frequency  Bin  Value 

freq2  (freqNum) 

Real 

Ending  Frequency  Bin  Value 

hsig 

Real 

Significant  Wave  Height 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 

ifreq 

Integer 

Number  of  Frequencies  in  Input  Spectrum 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

xffom  (dirNum) 

Real 

Direction  Array,  Direction  Wave  Energy 
Comes  From 

Output  Variables: 

esowm  (dirNum, freqNum) 

Real 

Directional  Wave  Spectrum 

Local  Variables: 

ang 

Real 

Temporary  Wave  .Angle 

b 

Real 

Constant  =  0.74 

const 

Real 

Variable  in  Pierson-Moskowitz  Equation 

e 

Real 

Variable  in  Pierson-Moskowitz  Equation 

enew 

Real 

Variable  in  Pierson-Moskowitz  Equation 

gu 

Real 

Variable  in  Pierson-Moskowitz  Equation 

hs 

Real 

Set  to  Significant  Wave  Height 

hsl 

Real 

Set  to  Significant  Wave  Height 

idir 

Integer 

Direction  Loop  Counter 

iffq 

Integer 

Frequency  Loop  Counter 

ipm 

Integer 

Set  to  1 

ratio 

Real 

Set  to  1.0 

sprd 

Real 

Directional  Spreading  Factor 

suml 

Real 

Temporary  Wave  Energy  Variable 

sum2 

Real 

Temporary  Wave  Energy  Variable 

temp 

Real 

Variable  in  Pierson-Moskowitz  Equation 

theta 

Real 

Wave  Angle 

vail 

Real 

Variable  in  Pierson-Moskowitz  Equation 

val2 

Real 

Variable  in  Pierson-Moskowitz  Equation 

wl 

Real 

Wave  Frequency  at  Beginning  of  Bin 

w2 

Real 

Wave  Frequency  at  End  of  Bin 

Subroutines  Called  from  SEAFIT  ( ): 

None. 
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SEAFIT  ( )  Called  from  Subroutines: 

WAVEFIT 

Figure  57.  Subroutine  SEAFIT  Flowchart 
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5.57  Subroutine  SETUP 


Subroutine  Call: 

SETUP  (pkfreq,  dl,  d2,  hrmsl,  hrms2,  etal,  eta2,  kinitl) 

Summary: 

Subroutine  SETUP  calculates  the  change  in  the  nearshore  mean  water  level  caused  by  the 
onshore  flux  of  momentum  or  the  shore-directed  Radiation  Stress.  The  presence  of  waves  causes  a 
change  in  the  total  water  depth,  which  is  defined  by  the  still  water  level  plus  the  wave-induced  set-up. 

Input  Variables: 


dl 

Real 

Corresponding  Depth 

d2 

Real 

Next  Corresponding  Depth 

etal 

Real 

Wave  Induced  Setup  at  Present  Location 

hrmsl 

Real 

Root  Mean  Square  Wave  Height 

hrms2 

Real 

Wave  Height  at  next  Onshore  Grid  Location 

kinitl 

Real 

Wave  Number 

pkfreq 

Real 

Peak  Frequency  at  the  Center  of  the 
Frequency  Band 

Output  Variables: 

eta2 

Real 

Wave  Induced  Setup  at  New  Location 

Local  Variables: 

avgdepth 

Real 

Averaged  Depth 

convrg 

Logical 

Set  to  False 

el 

Real 

Total  Average  Energy  for  Offshore  Wave 

e2 

Real 

Total  Average  Energy  for  Wave  Shoaled  and 
Refracted  Toward  the  Shore 

enl 

Real 

Linear  Wave  Theory  Ratio  of  Group 

Velocity  to  Wave  Celerity 

en2 

Real 

Linear  Wave  Theory  Ratio  of  Group 

Velocity  to  Wave  Celerity 

etanew 

Real 

Wave  Induced  Setup  Estimated  at 

New  Location 

i 

Integer 

Counter 

kl 

Real 

First  Wave  Number  Estimate 

k2 

Real 

Second  Wave  Number  Estimate 

percentdiff 

Real 

Convergence  Check 
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sxxl  Real  Cross-Shore  Directed  Radiation  Stress 

SXx2  Real  Cross-Shore  Directed  Radiation  Stress 

tol  Real  Convergence  Tolerance 

Subroutines  Called  from  SETUP  ( ): 

WAVENUM 

SETUP  ( )  Called  from  Subroutines: 

MAIN  WAV 
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Figure  58.  Subroutine  SETUP  Flowchart 
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5.58  Subroutine  SHORTOUT 


Subroutine  Call: 

SHORTOUT  (wdir,  wspd,  j,  iimax,  dxy,  xtemp,  suml,  k,  hlmax,  h2max,  per.  pet,  thetal, 
vmax,  vmin,  width,  igamma,  bl,  rk,  htemp,  widii,  j  gamma,  alfa,  bravo,  chrlie,  echo,  foxtrt,  golfl, 
golf2,  ihtll,  ihtl2) 

Summary: 

Subroutine  SHORTOUT  defines  the  forecasting  output  variables. 

Input  Variables: 


bl  (points) 

Real 

dxy  (points) 

Real 

hlmax 

Real 

h2max 

Real 

htemp  (points) 

Real 

igamma 

Integer 

iimax 

Integer 

j 

Integer 

k 

Integer 

pct(4) 

Real 

per 

Real 

rk  (points,  4) 

Real 

suml 

Real 

thetal 

Real 

vmax 

Real 

vmin 

Real 

wdir 

Real 

widii 

Integer 

width 

Real 

wspd 

Real 

xtemp  (points) 

Real 

Bottom  Slope  Array 
Corresponding  Depths  with  Tide 
Largest  Significant  Wave  Height  in  the 
Surf  Zone 

Largest  Maximum  Wave  Height  in  the 
Surf  Zone 

Temporary  Variable  for  Significant  Wave 
Height  Values 

Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Tow  ard  Beach 
Number  of  Calculation  Locations 
Pre-tidal  Depth  or  Still  Water  Level 
Temporary  Variable  for  Significant 
Wave  Height 
Percentage  Breaker  Array 
Peak  Period  of  Directional  Wave  Spectrum 
Matrix  of  Percentage  Breakers  and 
Type  of  Breakers 

Sum  of  Wave  Length  in  the  Surf  Zone 
Wave  Angle  at  Input  Starting  Depth 
Maximum  Positive  Longshore 
Current  Velocity 
Maximum  Negative  Longshore 
Current  Velocity 

Input  Wind  Direction  -  Compass  Heading 

Wind  is  Blowing  From 

Surf  Zone  Width  Array  Index 

Surf  Zone  Width 

Input  Wind  Speed 

Temporary  Variable  for  Cross-Shore  Values 


Output  Variables: 
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alfa 

Real 

Significant  Breaker  Height 

bravo 

Real 

Maximum  Breaker  Height 

chrlie 

Real 

Dominant  Breaker  Period 

echo 

Real 

Breaker  Angle 

foxtrt 

Real 

Longshore  Current  Speed  and  Direction 

golfl 

Real 

Number  of  Surf  Lines 

gol£2 

Real 

Surf  Zone  Width 

ihtll 

Real 

Wind  Speed 

ihtl2 

Real 

Wind  Direction 

jgamma 

Integer 

Temporary  Variable  Set  to 

Beach  Orientation 

Local  Variables: 

il 

Integer 

Temporary  Array 

i2 

Integer 

Temporary  Array 

tempi 

Real 

Temporary  Variable  for  Longshore  Current 
Maximum  Calculation 

temp2 

Real 

Temporary  Variable  for  Longshore  Current 
Minimum  Calculation 

xlen 

Real 

Average  Wave  Length  in  Surf  Zone 

Subroutines  Called  from  SHORTOUT( ): 

NEW_BRK 

SHORTOUT  ( )  Called  from  Subroutines: 

CALC SURF 
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Figure  59.  Subroutine  SHORTOUT  Flowchart 
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5.59  Subroutine  SLF  STRT 


Subroutine  Call: 

SLF_STRT  (self_st,  xk,  theta,  xdelt_gr,  hrms,  per.  fqz,  fqd,  Cg,  dxy,  nnn,  b, 
j  ii,  10,  thetaO,  surf) 

Summary: 

Subroutine  SLF_STRT  shoals  and  refracts  waves  from  the  farthest  offshore  point  to  the 
shoreward  point  where  the  percentage  of  breaking  exceeds  the  surf  zone  criteria  of  five  percent  (5%). 
If  the  five  percent  (5%)  threshold  is  not  exceeded,  execution  halts. 

Input  Variables: 


b 

Real 

Empirical  Factor  in  Wave  Breaking  Model 

Cg 

Real 

Wave  Group  Velocity 

dxy  (points) 

Real 

Corresponding  Depths  with  Tide 

fqd 

Real 

Peak  Frequency  at  the  Center  of  the 

Frequency  Band 

fqz 

Real 

Zero  Crossing  Frequency 

hrms 

Real 

Root  Mean  Square  Wave  Height 

nnn 

Integer 

Number  of  Points  in  Input  Depth  Array 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

self_st 

Char*  1 

Self  Staring  Option  (Yes  or  No) 

theta 

Real 

Radiation  Stress  Angle 

W) 

1 

o 

X 

Real 

Self-Adjusting  Cross-Shore  Grid  Step 

xk 

Real 

Wave  Number 

Output  Variables: 

Cg 

Real 

Wave  Group  Velocity 

hrms 

Real 

Root  Mean  Square  Wave  Height 

jji 

Integer 

Index  where  Wave  Probabilities 

Exceed  Threshold 

10 

Real 

Wave  Length  Offshore  Location 

surf 

Logical 

Index  Where  Percentage  of  Breakers  Is  Exceeded 
-  Start  of  Surf  Zone 

thetaO 

Real 

Wave  Angle  at  Grid  Offshore  Location 

xk 

Real 

Wave  Number 

Page  148 


Local  Variables: 


beta 

Real 

Bottom  Slope 

cg2 

Real 

Group  Velocity 

convg 

Real 

Convergence  Flag  (True  or  False) 

dp 

Real 

Offshore  Water  Depth 

eb 

Real 

Dissipation  Term 

hrms2 

Real 

Root  Mean  Square  Wave  Height 

ii 

Integer 

Array  Index 

1 

Real 

Wave  Length 

P(4) 

Real 

Breaker  Percentage  Array 

rhs 

Real 

Right  Hand  Side  of  Energy  Equation 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

rstart 

Real 

Percent  Breaking  Wave  Criteria 

xkO 

Real 

Offshore  Wave  Number 

Subroutines  Called  from  SLF_STRT  (  ): 

BALANCEQ 

GET_RHS 

PERCENT 

PT2 

SLF  STRT  ( )  Called  from  Subroutines: 

MAINJWAV 
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Figure  60.  Subroutine  SLF  STRT  Flowchart 
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5.60  Subroutine  SPLINE 


Subroutine  Call: 


SPLINE  (xi,  c,  n) 


Summary: 

Subroutine  SPLINE  calculates  the  coefficients  of  the  cubic  polynomial  that  fits  through  a 


specific  set  of  x  and  y  coordinates. 

Input  Variables: 


xi  (dirNum)  Real 

Output  Variables: 

c  (4, dirNum)  Real 

n  Integer 

Local  Variables: 

d  (dirNum)  Real 

dngl  (dirNum)  Real 

dvdl  Real 

dvd3  Real 

dx  Real 

gg  Real 

I  Integer 

m  Integer 


Subroutines  Called  from  SPLINE  ( ): 


Array  of  X-Coordinates 


Cubic  Polynomial  Coefficients 
Number  of  X-Coordinates 


Difference  between  Adjacent  X-Coordinates 
Slope  between  (2)  Adjacent  X-Coordinates 
Temporary  Variable 
Temporary  Variable 
Difference  between  (2)  points 
Temporary  Value  (delta  x  /  slope) 

Loop  Counter 
Loop  Counter 

None. 


SPLINE(  )  Called  from  Subroutines: 

GENRLSPL 
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Figure  61.  Subroutine  SPLINE  Flowchart 
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Number  ofX- 
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i 
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End  of  Loop 


\ 

RETURN  ] 
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5.61  Subroutine  SRFSETUP 


Subroutine  Call: 

SRFSETUP  (filein,  fileout,  fracname,  lndname,  depname,  iyear,  imonth,  iday,  ihour,  imin, 
gamma2,  ydepth,  slope,  ydetail,  xdelt,  dstart,  yrefrac,  ystr.  self_st,  hsea,  psea,  dsea,  hswell,  pswell, 
dswell,  wspd,  wdir,  tide,  gt_frg,  spefile,  file_dat,  file  tmp.  spedepth,  file_spc) 

Summary: 

Subroutine  SRFSETUP  opens  input  and  output  files.  Input  variables  are  initialized  using  data 
from  user-constructed  input  file.  The  format  of  the  input  file  is  outlined  in  Section  6.0. 

Input  Variables:  None. 


Output  Variables: 


depname 

Char*40 

Depth  Profile  File  Name 

dsea 

Real 

Input  Direction  for  Sea  Contribution 

dstart 

Real 

Input  Starting  Depth 

dswell 

Real 

Input  Swell  Direction  for  Internally 

Generated  Spectrum 

filein 

Char*40 

Input  File  Name 

file_out 

Char*  40 

Output  File  Name 

file_dat 

Char*  40 

Output  File  Name 

file_spc 

Char*  40 

Shallow  Water  Wave  Spectrum  File  Name 

file_tmp 

Char*  40 

Output  File  Name 

fracname 

Char*  40 

Wave  Refraction  File  Name 

gamma2 

Real 

Beach  Orientation,  Compass  Heading 

Directly  Toward  Beach 

gt_frg 

Integer 

Spectrum  Type 

hsea 

Real 

Input  Significant  Wave  Height  for  Sea 
Contribution  to  Pierson  Moskowitz  Spectrum 

hswell 

Real 

Input  Significant  Wave  Height  for  Internally 
Generated  Spectrum 

iday 

Integer 

Input  Day 

ihour 

Integer 

Input  Hour 

imin 

Integer 

Input  Minute 

imonth 

Integer 

Input  Month 

iyear 

Integer 

Input  Year 

lndname 

Char*  40 

Input  Landing  Zone  Name 

psea 

Real 

Input  Wave  Period  for  Sea  Contribution  to 
Internally  Generated  Spectrum 

pswell 

Real 

Input  Swell  Period  for  Internally 

Generated  Spectrum 
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selfst 

Char*  1 

slope 

Real 

spedepth 

Real 

spefile 

Char*40 

tide 

Real 

wdir 

Real 

wspd 

Real 

xdelt 

Real 

ydepth 

Char*  1 

ydetail 

Char*  1 

yrefrac 

Char*  1 

ystr 

Char*  1 

Local  Variables: 

duml 

Char*  80 

fend 

Integer 

filedat 

Char*20 

i 

Integer 

iopen 

Integer 

j 

Integer 

Self  Start  Flag  (Yes  or  No) 

Bottom  Slope 

Depth  at  Offshore  Wave  Spectrum 
Selected  Wave  Spectrum  File  Name 
Input  Tide  Level 

Input  Wind  Direction,  Compass  Heading 

Wind  Blows  From 

Input  Wind  Speed 

Surf  Zone  Output  Interval 

Input  Depth  Profile  Used?  (Yes  or  No) 

Detailed  Output?  (Yes  or  No) 

Is  Refraction  Considered  in  Analysis? 

(Yes  or  No) 

Is  Straight  Coast  Refraction  Used?  (Yes  or  No) 


Title  Line 

File  Name  Prefix  Used  for  Building 
File  Names 

Additional  Output  File  Name 
Loop  Counter 
I/O  Status  Number 
Loop  Counter 


Subroutines  Called  from  SRFSETUP  ( ): 

ABORT 

SRFSETUP  ( )  Called  from  Subroutines: 

SURF 
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Figure  62.  Subroutine  SRFSETUP  Flowchart 


I 


I 
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5.62  Subroutine  STRFRAC 


Subroutine  Call: 

STRFRAC  (dstart,  ifreq,  freq,  igamma,  idirec,  xfrom,  xcoeff,  xtheta,  wavedep) 

Summary: 

Subroutine  STRFRAC  calculates  wave  angle  refraction  coefficients  and  combined  shoaling 
and  refraction  coefficients  to  propagate  wave  energy  into  shallow  water. 


Input  Variables: 


dstart 

Real 

Input  Starting  Depth 

freq  (freqNum) 

Real 

Input  Wave  Spectrum  Center  Frequencies 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 

ifreq 

Integer 

Number  of  Frequency  Bins  in 

Input  Spectrum 

igamma 

Integer 

Beach  Orientation  Rotated  90  Degrees  from 
Original  Heading  Toward  Beach 

xfrom  (dirNum) 

Real 

Direction  Array 

Output  Variables: 

xcoeff  (dirNum, freqNum) 

Real 

Wave  Height  Refraction  Coefficients 

xtheta  (dirNum, freqNum) 

Real 

Wave  Angle  Refraction  Coefficients 

Local  Variables: 

argl 

Real 

Shallow  Water  Angle  (1)  -  Temporary 

direc 

Real 

Temporary  Direction  Angle 

frd 

Real 

Wave  Frequency 

idir 

Integer 

Direction  Loop  Counter 

ifrq 

Integer 

Frequency  Loop  Counter 

m 

Integer 

Temporary  Wave  Angle 

noprint 

Real 

Wave  Component  Direction 

shoal 

Real 

Temporary  Shoaling  Coefficient 

shoal2 

Real 

Temporary  Shoaling  Coefficient  at  Input 
Starting  Depth 

thetad 

Real 

Temporary  Wave  Angle  Variable 

thetas2 

Real 

Temporary  Wave  Angle  Variable 

xkd 

Real 

Temporary  Wave  Number  Variable 

xk2 

Real 

Temporary  Wave  Number  Variable 
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xks2 


Real 


Temporary  Wave  Number  at  Input 
Starting  Depth 

xksd2  Real  Wave  Number  at  Input  Starting  Depth 

Subroutines  Called  from  STRFRAC  ( ): 

WAVENUM 

STRFRAC  (  )  Called  from  Subroutines: 

SURF 
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Figure  63.  Subroutine  STRFRAC  Flowchart 
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5.63  Subroutine  SUMMARY 


Subroutine  Call: 

SUMMARY  (dstart,  tide,  wspd,  wdir,  xdelt,  yrefrac,  ystr,  depname,  file_out,  fracname, 
lndname,  ydepth,  ydetail,  gamma2,  slope,  hsea,  psea,  dsea,  hswell,  pswell,  dswell,  spectra,  spefile) 

Summary: 

Subroutine  SUMMARY  summarizes  the  input  information  read  to  the  output  file  for 
documentation  and  forecaster  verification. 

Input  Variables: 


depname 

Char*  40 

dsea 

Real 

dstart 

Real 

dswell 

Real 

fileout 

Char*  40 

fracname 

Char*  40 

gamma2 

Real 

hsea 

Real 

hswell 

Real 

lndname 

Char*  40 

psea 

Real 

pswell 

Real 

slope 

Real 

spectra 

Logical 

spefile 

Char*40 

tide 

Real 

wdir 

Real 

wspd 

Real 

xdelt 

Real 

ydepth 

Char*l 

ydetail 

Char*  1 

yrefrac 

Char*  1 

ystr 

Char*l 

Depth  Profile  File  Name 

Input  Direction  for  Sea  Contribution 

Input  Starting  Depth 

Input  Swell  Direction  for  Internally 

Generated  Spectrum 

Output  File  Name  *.out 

Wave  Refraction  File  Name 

Beach  Orientation,  Compass  Heading  Directly 

Toward  Beach 

Input  Significant  Wave  Height  for  Sea 

Contribution  to  Internally  Generated  Spectrum 

Input  Significant  Wave  Height  to  Internally 

Generated  Spectrum 

Input  Landing  Zone  Name 

Input  Wave  Period  for  Sea  Contribution  to 

Internally  Generated  Spectrum 

Input  Swell  Period  for  Internally 

Generated  Spectrum 

Bottom  Slope  for  a  Constructed  Depth  Profile 
Does  Input  Spectrum  Exist?  (True  or  False) 
Selected  Wave  Spectrum  File  Name 
Input  Tide  Level 

Input  Wind  Direction  Compass  Heading  Wind 

Blows  From 

Input  Wind  Speed 

Surf  Zone  Output  Interval 

Input  Depth  Profile  Used?  (Yes  or  No) 

Detailed  Output?  (Yes  or  No) 

Is  Refraction  Considered  in  Analysis? 

(Yes  or  No) 

Is  Straight  Coast  Refraction  Used?  (Yes  or  No) 
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Output  Variables: 


None. 


Local  Variables: 

sediment  Char*40 

Subroutines  Called  from  SUMMARY  ( ): 
SUMMARY  ( )  Called  from  Subroutines: 

SURF 


Sediment  Type 
None. 


Page  160 


Figure  64.  Subroutine  SUMMARY  Flowchart 
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5.64  Subroutine  SURFCAST 


Subroutine  Call: 

SURFCAST  (pet,  depname,  lndname,  slope,  ydepth.  alfa,  bravo,  chrlie,  echo,  foxtrt,  golfl, 
golf2,  ihtll,  ihtl2) 

Summary: 

Subroutine  SURFCAST  reads  input  variables  and  provides  a  short  format  summary  of  Navy 
specified  parameters.  The  subroutine  also  examines  longshore  current  direction  and  selects  the 
dominant  breaker  type. 


Input  Variables: 


alfa 

Real 

Significant  Breaker  Height 

bravo 

Real 

Maximum  Breaker  Height 

chrlie 

Real 

Dominant  Breaker  Period 

depname 

Char*40 

Depth  Profile  File  Name 

echo 

Real 

Breaker  Angle 

foxtrt 

Real 

Longshore  Current  Speed  and  Direction 

golfl 

Real 

Number  of  Surf  Lines 

golf2 

Real 

Surf  Zone  Width 

ihtll 

Real 

Wind  Speed  Coded  Surf  Forecast  Value 

ihtl2 

Real 

Wind  Direction 

lndname 

Char*40 

Input  Landing  Zone  Name 

pet  (4) 

Real 

Percent  of  Different  Breaker  Types: 

pet  (1)  =  Spilling 

pet  (2)  =  Plunging 

pet  (3)  =  Surging 

pet  (4)  =  Total 

slope 

Real 

Bottom  Slope 

ydepth 

Char*l 

Input  Depth  Profile  Used?  (Yes  or  No) 

Output  Variables: 

None. 

Local  Variables: 

foxtmp 

Real 

Longshore  Current  Where  the  Sign  Indicates  the 
Direction 

i 

Integer 

Loop  Counter  Variable 

jdelt 

Integer 

Difference  If  Any  Between  100%  and 

Sum  ofjp  (4) 

jP  (4) 

Integer 

Temporary  Variable  Same  as  pct(4)  Array 
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jsum 

Integer 

Check  for  Percentages  Adding  to  100% 

maxp 

Integer 

Indicates  Dominant  Breaker  Type 

xmax 

Real 

Temporary  Variable  Used  in  Dominant  Breaker 
Type  Examination 

Subroutines 

Called  from  SURFCAST  ( ): 

None. 

SURFCAST 

( )  Called  from  Subroutines: 

SURF 
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Figure  65.  Subroutine  SURFCAST  Flowchart 


[  START 


I 


Determine  Wave 
Occurrence 
Percentages 


I 


Output  Breaker 
Types 


I 


Verify  Wav  e 
Percentages 
=  100% 


Output  Littoral 
Current  Direction 


I 


Page  164 


5.65  Subroutine  SWLFIT 


Subroutine  Call: 

SWLFIT  (hsig,  per,  dir,  dangle,  ifreq,  idirec,  period,  esowm) 

Summary: 

Subroutine  SWLFIT  superimposes  remotely  generated  swell  wave  energy  onto  the  existing 
directional  wave  spectrum.  The  existing  wave  spectrum  may  be  zero  or  it  may  contain  locally 
generated  sea  waves  already  added  by  the  subroutine  SEAFIT. 


Input  Variables: 


dangle 

Real 

Angle  between  Directional  Bins 

dir 

Real 

Input  Swell  Direction  for  Internally 

hsig 

Real 

Generated  Spectrum 

Significant  Wave  Height 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 

ifreq 

Integer 

Number  of  Frequency  Bins  in 

per 

Real 

Input  Spectrum 

Peak  Period  of  Directional  Wave  Spectrum 

period  (freqNum) 

Real 

Period  Array  (1  /  Frequency) 

Output  Variables: 

esowm  (dirNum, freqNum) 

Real 

Directional  Wave  Spectrum 

Local  Variables: 

dl 

Real 

Temporary  Variable  for  Distributing 

d2 

Real 

Wave  Energy 

Temporary  Variable  for  Distributing 

d3 

Real 

Wave  Energy 

Temporary  Variable  for  Distributing 

delt 

Real 

Wave  Energy 

Temporary  Variable  for  Distributing 

diff 

Real 

Wave  Energy 

Difference  between  Maximum  Wave 

dmin 

Real 

Period  and  Array  Value  of  Wave  Period 

Set  to  1000.0 

energy 

Real 

Swell  Energy 

ifrq 

Integer 

Frequency  Loop  Counter 

jdir 

Integer 

Swell  Direction 
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jdirl 

Integer 

jdir3 

Integer 

jfreq 

Integer 

xdir 

Real 

Subroutines  Called  from  SWLFIT  ( ): 
SWLFIT  ( )  Called  from  Subroutines: 

WAVEFIT 


Figure  66.  Subroutine  SWLFIT  Flowchart 


Direction  Bin  Index  Number 
Direction  Bin  Index  Number 
Directional  Wave  Spectrum  Wave  Number 
Wave  Direction 

None. 
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5.66  Subroutine  WAVEFIT 


Subroutine  Call: 

WAVEFIT  (ifreq,  idirec,  dangle,  hsea,  psea,  dsea,  hswell,  pswell,  dswell,  freql,  freq2,  xfrom, 
period,  esowm,  ehsig) 

Summary: 

Subroutine  WAVEFIT  initializes  the  internally  generated  directional  wave  spectrum  to  zero 
and  calls  subroutines  SEAFIT  and  SWLFIT  to  fill  the  matrix. 


Input  Variables: 


dangle 

Real 

Angle  Between  Directional  Bins 

dsea 

Real 

Input  Direction  for  Sea  Contribution 
to  Internally  Generated  Wave  Spectrum 

dswell 

Real 

Input  Swell  Direction  for  Internally 
Generated  Spectrum 

freql  (freqNum) 

Real 

Beginning  Frequency  Bin  Value 

freq2  (freqNum) 

Real 

Ending  Frequency  Bin  Value 

hsea 

Real 

Input  Significant  Wave  Height  for  Sea 
Contribution  to  Internally  Generated 

Wave  Spectrum 

hswell 

Real 

Input  Significant  Wave  Height  to 

Internally  Generated  Spectrum 

idirec 

Integer 

Number  of  Direction  Bins  in  Input  Spectrum 

ifreq 

Integer 

Number  of  Frequencies  in  Input  Spectrum 

period  (freqNum) 

Real 

Period  Array  (1/Frequency) 

psea 

Real 

Input  Wave  Period  for  Sea  Contribution 

pswell 

Real 

Input  Swell  Period  for  Internally 

Generated  Spectrum 

xfrom  (dirNum) 

Real 

Direction  Array,  Direction  Wave  Energy 
Comes  From 

Output  Variables: 

ehsig 

Real 

Significant  Wave  Height  from 

Directional  Spectrum 

esowm  (dirNum,  freqNum) 

Real 

Directional  Wave  Spectrum 
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Local  Variables: 

idir  Integer  Direction  Loop  Counter 

ifrq  Integer  Frequency  Loop  Counter 

Subroutines  Called  from  WAVEFIT  ( ): 

GTSIGH 

SEAFIT 

SWLFIT 

WAVEFIT  ( )  Called  from  Subroutines: 

GENSPEC 

Figure  67.  Subroutine  WAVEFIT  Flowchart 
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5.67  Subroutine  WAVENUM 


Subroutine  Cali: 

WAVENUM  (fq,  dp,  xk) 

Summary: 

The  wave  dispersion  equation  is  solved  for  the  wave  number  through  numerical  iteration.  A 
relative  change  of  less  than  .0005  is  required  and  the  maximum  number  of  iterations  is  150.  If 
convergence  is  not  obtained  within  150  iterations,  a  shallow  water  approximation  is  employed. 


Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

fq 

Real 

Frequency  either  Wave  or  Peak 

Output  Variables: 

xk 

Real 

Wave  Number 

Local  Variables: 

const 

Real 

Shallow  Water  Criteria  Constant 

diff 

Real 

Percent  Difference  between  Wave 
Number  Estimates 

est 

Real 

Estimate  of  Wave  Number 

I 

Integer 

Loop  Counter 

it 

Integer 

Loop  Limit  -  Set  to  150 

Subroutines  Called  from  WAVENUM  ( ): 

None. 

WAVENUM  ( )  Called  from  Subroutines: 

INITLIZE 

PT2 

RADST1 

RADST2 

SETUP 

STRFRAC 
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Figure  68.  Subroutine  WAVENUM  Flowchart 
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5.68  Subroutine  WEIGHTFN 


Subroutine  Call: 

WEIGHTFN  (dp,  hrms,  h,  w_h) 

Summary: 

Subroutine  WEIGHTFN  calculates  the  weighting  function  used  to  describe  the  distribution 
of  breaking  waves  across  the  surf  zone. 

Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

h 

Real 

Wave  Height 

hrms 

Real 

Root  Mean  Square  Wave  Height 

Output  Variables: 

w_h 

Real 

Output  Weighting  Function 

Local  Variables: 

m 

Real 

Multiplier 

temp 

Real 

Weighting  Function 

tol 

Real 

Set  to  -700.00 

Subroutines  Called  from  WEIGHTFN  (  ): 

GETM 

WEIGHTFN  ( )  Called  from  Subroutines: 

F2 
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Figure  69.  Subroutine  WEIGHTFN  Flowchart 
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5.69  Subroutine  ZONE1 


Subroutine  Call: 

ZONE1  (jji,  iimax,  dxy,  xtemp,  htemp,  ptemp,  xktemp,  v,  distmax,  vmax,  vmin,  suml, 
width,  j,  k,  hlmax,  h2max,  widii) 

Summary: 

Subroutine  ZONE1  calculates  the  preliminary  surf  forecast  values  and  surf  zone  parameters. 


Input  Variables: 

distmax 

Real 

Farthest  Distance  Offshore 

dxy  (points) 

Real 

Pre-Tidal  Depth  or  Still  Water  Level 

htemp  (points) 

Real 

Temporary  Variable  for  Significant  Wave 
Height  Values 

iimax 

Integer 

Number  of  Calculation  Locations 

jji 

Integer 

Index  where  Wave  Probabilities 

Exceed  Threshold 

ptemp  (points) 

Real 

Percentage  of  Breaking  Waves  and 

Breaker  Types 

v  (points) 

Real 

Longshore  Current 

xktemp  (points) 

Real 

Temporary  Variable  for  Wave  Number 

xtemp  (points) 

Real 

Temporary  Variable  for  Cross-Shore  Values 

Output  Variables: 

hlmax 

Real 

Maximum  Significant  Wave  Height 

h2max 

Real 

Maximum  Wave  Height 

j 

Integer 

Array  Index  Where  Maximum  Significant  Wave 
Height  Occurs 

k 

Integer 

Temporary  Variable  Number  of  Points  in  Cross- 
Shore  Transect 

suml 

Real 

Summation  of  Wave  Lengths  Across  the 

Surf  Zone 

vmax 

Real 

Maximum  Positive  Longshore  Current  Velocity 

vmin 

Real 

Maximum  Negative  Longshore 

Current  Velocity 

widii 

Integer 

Array  Index  for  X- value  at  Surf 

Zone  Boundary 

width 

Real 

Surf  Zone  Width 
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Local  Variables: 


dpi 

Real 

hdep 

Real 

hmax 

Real 

houtl 

Real 

hrmsl 

Real 

ii 

Integer 

vingl 

Real 

wlen 

Real 

xoffl 

Real 

Offshore  Depth  in  Feet 
Limiting  Breaking  Depth 
Temporary  Variable  for 
Maximum  Wave  Height 
Temporary  Variable  for  Significant 
Wave  Height 

Root  Mean  Square  Wave  Height 
Loop  Index 

Longshore  Current  Velocity  in  Knots 
Wave  Length 
Distance  Offshore 


Subroutines  Called  from  ZONE1  ( ):  None. 

ZONE1  (  )  Called  from  Subroutines: 


CALCSURF 


Figure  70.  Subroutine  ZONE1  Flowchart 


5.70  Function  CUBPOLY 


Function  Call: 


CUBPOLY  (xavg,  xi,  c,  n) 

Summary: 


Function  CUBPOLY  evaluates  the  cubic  polynomial  that  was  previously  fit  through  a  defined 


set  of  x  and  y  coordinates.  The  evaluated  cubic 
an  input  x  value. 

Input  Variables: 


c  (4,dirNum)  Real 

n  Integer 

xavg  Real 

xi  (dirNum)  Real 

Output  Variables: 

CUBPOLY  Real 

Local  Variables: 

dx  Real 

I  Integer 

j  Integer 


polynomial  function  interpolates  a  new  y  value  for 


Cubic  Polynomial  Coefficient 
Number  of  X-Coordinates 
Interpolated  Coordinate 
Array  of  X-Coordinate 


Value  at  the  Interpolated  Coordinate 


Temporary  Variable 
Loop  Counter 
Loop  Counter 


Subroutines  Called  from  CUBPOLY  (  ):  None. 

CUBPOLY  ( )  Called  from  Subroutines: 

GENRLSPL 
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Figure  71.  Function  CUBPOLY  Flowchart 
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5.71  Function  F2 


Function  Call: 

F2  (h,  hrms,  dp,  p  flag) 

Summary: 

Function  F2  evaluates  the  Rayleigh  probability  distribution  function  for  a  given  wave  height 
value,  for  a  selected  weighting  function. 

Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

h 

Real 

Wave  Height 

hrms 

Real 

Root  Mean  Square  Wave  Height 

pflag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

Output  Variables: 

£2 

Real 

Weighted  Rayleigh  Distribution 

Local  Variables: 

P_h 

Real 

Rayleigh  Probability  Distribution 

temp 

Real 

Exponent  Term  in  Rayleigh  Distribution 

tol 

Real 

Tolerance  Value  Set  to  -700.0 

w_h 

Real 

Weighting  Function 

Subroutines  Called  from  F2  ( ): 

WEIGHTFN 

F2  ( )  Called  from  Functions: 
INTEGRAT 
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Figure  72.  Function  F2  Flowchart 


I 
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5.72  Function  F3 


Function  Call: 

F3  (hrms,  theta,  Cg.  dp,  mean_freq,  xk,  roller) 

Summary: 

Function  F3  returns  values  for  the  LHS  of  the  energy  equation. 

Input  Variables: 


Cg 

Real 

Wave  Group  Velocity 

dp 

Real 

Offshore  Water  Depth 

hrms 

Real 

Root  Mean  Square  Wave  Height 

mean_freq 

Real 

Directional  Spectrum  Value 

roller 

Logical 

Roller  Option  Flag  (True  or  False) 

theta 

Real 

Wave  Angle 

xk 

Real 

Wave  Number 

Output  Variables: 

G 

Real 

Total  Energy 

Local  Variables: 

e_roller 

Real 

Roller  Contribution  to  the  Energy  Equation 

e_wave 

Real 

Wave  Contribution  to  the  Energy  Equation 

Subroutines  Called  from  F3  ( ): 

CALCROLL 

GETWAVE 

F3  ( )  Called  from  Subroutines: 

BALANCEQ 
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Figure  73.  Function  F3  Flowchart 
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5.73  Function  FCN1 


Function  Call: 

FCN1  (t,  h,  per,  1,  dp,  v,  u,  theta2 ) 

Summary: 

Function  FCN1  calculates  the  nonlinear  bottom  stress  term  at  a  particular  time.  This 
calculation  ultimately  provides  the  bottom  friction  for  the  longshore  current  calculation  after 
time-averaging  over  one  wave  period. 


Input  Variables: 


dp 

Real 

Offshore  Water  Depth 

h 

Real 

Wave  Height 

1 

Real 

Wave  Length 

per 

Real 

Peak  Period  of  Directional  Wave  Spectrum 

t 

Real 

Wave  Period 

theta2 

Real 

Wave  Angle 

u 

Real 

Mean  Cross-Shore  Current  Velocity 

V 

Real 

Longshore  Current  Velocity 

Output  Variables: 

fcnl 

Real 

Nonlinear  Bottom  Friction  at  a 

Specific  Time 

Local  Variables: 

d2 

Real 

Temporary  Variable  Used  in  Calculation 

uw 

Real 

Orbital  Velocity  at  Specific  Time 

um 

Real 

Orbital  Velocity 

w 

Real 

Angular  Frequency 

Subroutines  Called  from  FCN1  ( ): 

GET_UM 

FCN1  ()  Called  from  Subroutines: 

GETFCN 

Figure  74.  Function  FCN1  Flowchart 
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5.74  Function  INTEGRAT 


Function  Call: 

INTEGRAT  (xo,  xn,  hrms,  dp,  p_flag) 

Summary: 

Function  INTEGRAT  evaluates  an  integral  numerically  using  the  trapezoidal  rule.  Function 
{F2}  is  called  to  evaluate  the  integral  at  upper  and  lower  limits.  The  function  applies  the  trapezoidal 
integration  method  to  estimate  the  wave  height  at  a  particular  depth  from  a  weighted  distribution. 

Input  Variables: 


dp 

Real 

Farthest  Offshore  Water  Depth 

hrms 

Real 

Root  Mean  Square  Wave  Height 

pflag 

Logical 

Weighting  Factor  Flag  (True  or  False) 

xn 

Real 

Upper  Limit  of  Integration  =  5  *  hrms 

xo 

Real 

Lower  Limit  of  Integration  =  0.0 

Output  Variables: 

integrat 

Real 

Wave  Height  Distribution  Calculated  for  a 
Specific  Location 

Local  Variables: 

delt 

Real 

Step  Between  Intervals 

f_xn 

Real 

f(x)  Evaluated  at  Upper  Limit 

f  xo 

Real 

f(x)  Evaluated  at  Lower  Limit 

a 

Real 

Wave  Height  Distribution 

Weighting  Function 

i 

Integer 

Loop  Variable 

numit 

Integer 

Set  to  100  -  Number  of  Iterations  Examined  Over 
Integral 

sum 

Real 

Summary  Results  from  Function  F2 

xi 

Real 

Integration  Step  Location 

Functions  Called  from  INTEGRAT  ( ): 

F2 
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5.75  Include  File:  COMMON.INC 


Summary: 

The  include  file  COMMON.INC  contains  all  the  parameters  set  for  the  SURF  Model. 


Defined  Parameters: 

deal 

Real 

0.3048  -  Feet  to  Meters  Conversion 

degrad 

Real 

PI  / 180.0  -  Conversion  from 

Degrees  to  Radians 

dirNum 

Integer 

1 80  -  Array  Dimension  Used  for 

Direction  Arrays 

freqNum 

Integer 

50  -  Array  Dimension  Used  for 

Frequency  Arrays 

g 

Real 

9.8 

gamma 

Real 

0.42  -  Empirical  Wave  Height  Factor 

iunit 

Integer 

Output  File  Unit 

Pi 

Real 

3.14159265 

points 

Integer 

500  -  Array  Dimension  Used  for  all  Input 
Depth  Arrays 

raddeg 

Real 

1 80.0  /  pi  -  Conversion  from 

Radians  to  Degrees 

rho 

Integer 

1030  -  Water  Density 

rhoair 

Real 

1 .2  -  Air  Density 

sigma 

Real 

sigma_deg  *  degrad 

sigma-deg 

Real 

5.0  -  Angle  in  Degrees  between  Wave  and 
Roller  in  the  Thomton/Lippman  Model  (1996) 

tpi 

Real 

2  *  3.14159265 

zone_pct 

Real 

10%  Surf  Zone  Width  Percent  of  Breaking 
Waves 
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6.0 


NOTES 


6.1  SURF  3.1  Input  Files 


6.1.1  SURF  3.1  Input  File 


Line 

Description 

Type 

Units 

Range 

Line  1 

Input  File  Name 

Char*40 

Line  2 

Date  and  Time 

Char*  10 

(YYYYMMDDHH) 

Line  3 

Landing  Zone  Name 

Char*40 

Line  4 

Input  Depth  Profile  File  Name 

Char*40 

*  * 

Line  5 

Input  Wave  Spectrum  File  Name 

Char*40 

*  * 

Line  6 

Input  Wave  Refraction  File  Name 

:  Char*40 

*  * 

Line  7 

Compass  Heading  Toward  Beach 

Real 

Degrees 

0-359 

Line  8 

Slope/Sediment  Type 

Integer 

1-10 

1  =  Boulders  6  = 

Coarse  Sand 

2  =  Cobble  7  = 

Medium  Sand 

3  =  Pebbles  8  = 

Fine  Sand 

4  =  Granules  9  = 

Very  Fine  Sand 

5  =  Very  Coarse  Sand  10  = 

Silt 

Line  9 

Starting  Depth 

Real 

Feet 

>0 

Line  10 

Offshore  Wave  Spectrum  Depth 

Real 

Feet 

>0 

Line  1 1 

Sea  Wave  Height 

Real 

Feet 

>0 

Sea  Wave  Period 

Real 

Seconds 

1-30 

Sea  Wave  Direction 

Real 

Degrees 

0-359 

Swell  Wave  Height 

Real 

Feet 

>0 

Swell  Wave  Period 

Real 

Seconds 

1-30 

Swell  Wave  Direction 

Real 

Degrees 

0-359 

Line  12 

Wind  Speed 

Real 

Knots 

>0 

Wind  Direction 

Real 

Degrees 

0-359 

Tide  Elevation 

Real 

Feet 

+  or- 

Line  13 

Output  Data  Grid  Spacing 

Real 

Feet 

>0 

*  The  input  file  name  (line  1 )  must  always  be  included. 

**  If  any  of  the  above  input  data  is  not  included  or  not  available, 

insert  a  blank  or  a  blank  line  f 

character  and/or  numeric  data  to  maintain  a  consistent  format  in  the  input  file. 


***  The  above  format  is  for  the  default  model  setup,  for  more  detailed  information  read  the  advanced 
user  options  information  in  section  6.1.5. 
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6.1.2 

Line 

SURF  3.1  Input  Depth  Profile  File 

Description 

Type 

Range 

Line  1 

Title 

Char*80 

Line  2 

Units  for  Distance  Offshore 

1  -  Distances  in  Feet 

2  -  Distances  in  Meters 

3  -  Distances  in  Yards 

Integer 

1,2, or  3 

Line  3 

Units  for  Depth 

1  -  Depths  in  Feet 

2  -  Depths  in  Meters 

3  -  Depths  in  Fathoms 

Integer 

1,2, or  3 

Line  4 

-  EOF  Point  Number 

Integer 

1  -500 

Distance 

(+)  Positive  numbers  are  Offshore 
(-)  Negative  numbers  are  Onshore 

Real 

+  or  - 

Depth 

(+)  Positive  numbers  are  Depths 
(-)  Negativ  e  numbers  are  Elevations 

Real 

+  or  - 
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6.1.3  SURF  3.1  Wave  Refraction  File 


Line _ Description _ Type _ Units _ Range 


Line  1 

Longitude 

Real 

Degrees 

-180.0- +180. 

Line  2 

Latitude 

Real 

Degrees 

-90.0-  +90. 

Line  3 

Date  (YYYYDDMM) 

Real 

— 

Line  4 

Number  of  Angles 

Integer 

1  - 180 

Line  5 

Number  of  Rows 

Integer 

+  number 

Number  of  Columns 

Integer 

+  number 

Line  6 

Number  of  Frequency  Bins 

Integer 

1-50 

Line  7 

Initial  Direction 

Real 

Degrees 

0.  -  359. 

Line  8 

Initial  Frequency  Bin 

Real 

Degrees 

0.  -  359. 

Line  9 

Width  of  Direction  Bin 

Real 

Degrees 

2.  - 180. 

Line  10 

Direction  of  Waves 

1  -  Direction  waves  are  coming  from 

2  -  Direction  waves  are  going  to 

Integer 

1  or  2 

Angle  Refraction  Coefficients  -  This  section  is  repeated  for  each  Frequency  Bin 

Line 

Bin  Number 

Intege  r 

1-50 

Lower  Limit  of  Frequency  Bin 

Real 

Hertz 

>  =0. 

Center  of  Frequency  Bin 

Real 

Hertz 

>  =0. 

Upper  Limit  of  Frequency  Bin 

Real 

Hertz 

>  =0. 

Line 

Angle  Refraction  Coefficients  Real 

The  coefficients  are  in  the  format: 

(Number  of  Rows  by  Number  of  Columns). 

Degrees 

0.  -  359. 

All  rows  and  columns  must  contain  numbers;  pad  with  zeros,  if  necessary. 

End  of  Angle  Refraction  Coefficients  Section 

Line 

Header  1  for  Shoaling  Coefficients 

Char*  80 

Line 

Header  2  for  Shoaling  Coefficients 

Char*  80 

Line 

Header  3  for  Shoaling  Coefficients 

Char*  80 

Shoaling  Coefficients  -  This  section  is  repeated  for  each  Frequency  Bin 

Line 

Bin  Number 

Integer 

1-50 

Lower  Limit  of  Frequency  Bin 

Real 

Hertz 

>  =  0. 

Center  of  Frequency  Bin 

Real 

Hertz 

>  =  0. 

Upper  Limit  of  Frequency  Bin 

Real 

Hertz 

>  =  0. 

Line 

Shoaling  Coefficients 

Real 

(N/m)2 

The  coefficients  are  in  the  format: 

(Number  of  Rows  by  Number  of  Columns). 

All  rows  and  columns  must  contain  numbers;  pad  with  zeros,  if  necessary 


End  of  Shoaling  Coefficients  Section 

*  The  coefficients  in  this  file  must  be  defined  over  the  entire  0  to  360  degree  range.  A  partial  sector 
definition  (e.g.  0  to  180  degrees)  will  cause  errors.  If  the  input  data  is  not  available  over  the  entire 
range,  pad  the  direction  bins  with  zeros. 
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6.1.4  SURF  3.1  Spectrum  File 


Line 

Description 

Type 

Units 

Range 

Line  1 

Longitude 

Real 

Degrees 

-180.0 -+180.0 

Line  2 

Latitude 

Real 

Degrees 

-90.0  -  +90.0 

Line  3 

Date  -  (YYYYMMDD) 

Real 

— 

Line  4 

Number  of  Angles 

Integer 

1  - 180 

Line  5 

Number  of  Rows 

Integer 

+  number 

Number  of  Columns 

Integer 

+  number 

Line  6 

Number  of  Frequency  Bins 

Integer 

1-50 

Line  7 

Initial  Direction 

Real 

Degrees 

0.-359. 

Line  8 

Initial  Frequency  Bin 

Real 

Hertz 

>  =  0. 

Line  9 

Width  of  Direction  Bin 

Real 

Degrees 

2.-180. 

Line  10 

Direction  of  Waves 

1  -  Direction  waves  are  coming  from 

2  -  Direction  waves  are  going  to 

Integer 

1  or  2 

Directional  Wave  Spectrum  -  This  section  is  repeated  for  each  Frequency  Bin 

Line 

Bin  Number 

Integer 

1  -50 

Lower  Limit  of  Frequency  Bin 

Real 

Hertz 

>  =  0 

Center  of  Frequency  Bin 

Real 

Hertz 

>  =  0 

Upper  Limit  of  Frequency  Bin 

Real 

Hertz 

>  =  0 

Line 

Directional  Wave  Spectrum  Real  > 

( 

The  Number  of  Angles  are  in  the  format: 

(Number  of  Rows  by  Number  of  Columns) 

m  ) 

Hz  *  Radians 

>  =  0 

All  rows  and  columns  must  contain  numbers;  pad 

fields  with  zeros,  if  necessary. 

End  of  Directional  Wave  Spectrum  Section 


*  The  coefficients  in  this  file  must  be  defined  over  the  entire  0  to  360  degree  range.  A  partial  sector 
definition  (e.g.  0  to  180  degrees)  will  cause  errors.  If  the  input  data  is  not  available  over  the  entire 
range,  pad  the  direction  bins  with  zeros. 


Page  190 


6.1.5  Advanced  SURF  3.1  Model  Options 

Several  run-time  model  options  included  in  Surf  3.1  are  transparent  to  the  user.  These  options 
are  reserved  for  the  advanced  or  expert  user  applying  the  model  to  unique  situations.  The  default  input 
settings  described  in  section  6.5.1  are  appropriate  for  most  model  runs.  However,  if  necessary  the  user 
can  control  the  wave  refraction  and  the  amount  of  output  data  including  the  production  of  an  additional 
file  with  a  shallow  water  wave  spectrum  after  transformation  due  to  shoaling  and  refraction.  These 
options  are  not  recommended  for  most  users. 

Wave  Refraction  Option 

The  default  wave  refraction  setting  uses  linear  wave  theory  and  Snell=s  Law  to  refract  waves 
with  a  straight  coast  assumption.  A  coast  is  assumed  straight  if  the  bottom  contours  are  straight  and 
generally  parallel  with  the  coastline.  Line  6  in  the  input  file  is  used  to  specify  an  externally  generated 
wave  refraction  file  that  includes  refraction  and  shoaling  coefficients.  Programs  such  as  REFDIF  and 
STWAVE  can  be  used  to  calculate  these  types  of  coefficients.  If  an  expert  user  wants  to  ignore  all 
refraction  effects  Line  6  must  contain  the  word  none  or  NONE.  This  option  is  not  recommended  for 
most  users. 

Self-Start  Option 

The  model  is  typically  configured  to  use  the  self-start  option.  This  option  expedites  model 
execution  by  shoaling  and  refracting  the  offshore  wave  spectrum  to  the  starting  depth  specified  in  Line 
9  of  the  input  file.  The  model  then  begins  stepwise  calculations  from  this  point  shoreward.  There  are 
two  advanced  user  options  associated  with  the  starting  depth.  These  options  are  selected  by  using  a 
negative  number  or  a  zero  in  Line  9  of  the  input  file. 

If  Line  9  of  the  input  file  contains  a  negative  number  the  self-start  option  will  not  be  used  and 
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the  data  written  to  the  output  file  will  begin  at  the  absolute  value  of  the  starting  depth  specified  in  Line 
9.  For  example,  if  Line  9  in  the  input  file  is  B15,  then  the  self-start  option  will  not  be  used  and  the 
columnar  data  in  the  output  file  will  begin  at  a  water  depth  of  1 5  feet.  If  Line  9  is  a  zero  the  self  start 
option  will  not  be  used  and  calculations  will  begin  at  the  farthest  point  offshore  as  defined  in  the  input 
depth  file  or  the  constructed  equilibrium  depth  profile. 

Wave  Spectrum  Depth  Option 

Line  10  of  the  input  file  is  used  to  specify  the  water  depth  at  the  input  directional  wave 
spectrum.  If  this  value  is  left  blank  or  is  defined  as  zero,  the  model  will  assume  that  the  wave 
spectrum  is  located  in  deep  water.  If  Line  10  of  the  input  file  is  a  negative  number,  an  additional 
output  file  is  created  with  the  shallow  water  directional  wave  spectrum  at  the  depth  specified  in  Line 
9,  the  starting  depth.  This  wave  spectrum  has  been  shoaled  and  refracted  to  the  starting  depth. 

The  format  of  this  ASCII  text  file  is  a  simple  matrix  of  rows  and  columns.  It  has  the  same 
name  as  the  output  file  except  that  the  file  extension  will  be  *.dws.  The  first  row  contains  the  center 
frequency  bin  definitions  and  the  first  column  defines  the  wave  direction  bins.  The  heart  of  the  matrix 
is  the  wave  energy  per  frequency  and  direction  with  the  units  m2/Hz*radians.  This  spectrum  has  the 
same  units  as  the  input  directional  wave  spectrum  so  that  direct  comparisons  can  be  made.  This  option 
is  available  for  users  interested  in  examining  the  transformation  of  the  directional  wave  spectrum  in 
shallow  water. 

Detailed  Output  Option 

The  final  advanced  user  option  controls  the  amount  of  data  in  the  output  file.  The  default 
option  will  create  an  output  file  with  the  detailed  output  of  columnar  data  of  many  wave  parameters 
across  the  surf  zone.  The  distance  between  each  of  these  points  is  defined  by  Line  13  in  the  input  file. 
If  Line  13  is  zero  or  a  negative  number,  only  the  summary  of  the  wave  parameters  in  the  coded  surf 
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forecast  will  be  reported  in  the  output  file  excluding  the  detailed  output. 
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6.2 


SURF  3.1  Output  Files 


6.2.1  SURF  3.1  Detailed  Output  File 

The  SURF  Detailed  Output  File  has  three  output  sections  delineated  by  asterisks.  The  first 
section  contains  the  input  parameters  and  several  variables  describing  the  directional  wave  spectrum. 
The  second  section  is  the  coded  surf  forecast  with  variables  specific  to  military  surf  observations.  The 
final  section  is  the  detailed  surf  output,  which  is  columnar  data  describing  cross-shore  distributions 
of  several  variables  including  wave  height,  water  depth,  wave  breaking,  and  longshore  current.  The 
filename  generated  is  “*.out”,  where  the  is  replaced  with  the  prefix  of  the  input  file  name. 


Line 

Description 

Type 

Units 

Line  1 

Surf  Header 

Character 

Line  2 

Blank  Line 

Character 

Line  3 

SURF  Model  Version 

Character 

Line  4 

Date  and  Time  of  Forecast 

Character 

Line  5 

Output  File  Name  Information 

Character 

Line  6 

Straight  Coast  Wave  Refraction  Option 

Character 

Line  7 

Landing  Zone  Name 

Character 

Line  8 

Sight  Line  Toward  Beach 

Real 

Degrees 

Line  9 

Interval 

Real 

Feet 

Line  10 

Starting  Depth 

Real 

Feet 

Line  1 1 

Depth  Profile  Name  or  Beach  Sediment  Type 

Character 

Line  12 

Spectrum  Usage  Text 

Or 

Character 

Sea  Wave  Height 

Real 

Feet 

Sea  Period 

Real 

Seconds 

Sea  Direction 

Real 

Degrees 

Line  13 

Spectrum  File  Name 

Or 

Character 

Swell  Wave  Height 

Real 

Feet 

Swell  Period 

Real 

Seconds 

Swell  Direction 

Real 

Degrees 

Line  14 

Wind  Speed 

Real 

Knots 

Line  15 

Wind  Direction 

Real 

Degrees 

Line  16 

Tide  Level 

Real 

Feet 

Line  17 

Blank  Line 

Character 
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Line 

Description 

TyPe _ 

Units 

Line  18 

Internal  Grid  Spacing 

Real 

Feet 

Line  19 

Significant  Wave  Height  from  Input  File 

Real 

Feet 

Line  20 

Significant  Wave  Height  from  Straight  Coast 

Real 

Feet 

Line  21 

Input  Spectrum  Type 

Integer 

Line  22 

Significant  Wave  Height  Offshore 

Real 

Feet 

Line  23 

Significant  Wave  Height 

Real 

Feet 

Line  24 

Peak  Frequency 

Real 

Hertz 

Line  25 

Zero-Crossing  Frequency 

Real 

Hertz 

Line  26 

Peak  Period 

Real 

Seconds 

Line  27 

Percentage  Breaking  Waves  at  Starting  Depth 

Real 

Percent 

Line  28 

Self  Starting  Option 

Character 

Line  29 

Blank  Line 

Character 

Line  30 

Text  Heading  -  Surf  Forecast 

Character 

Line  3 1 

Significant  Breaker  Height 

Real 

Feet 

Line  32 

Maximum  Breaker  Height 

Real 

Feet 

Line  33 

Dominant  Breaker  Period 

Real 

Seconds 

Line  34 

Dominant  Breaker  Type 

Character 

Line  35 

Breaker  Percentages 

Character 

Percent 

Line  36 

Breaker  Angle 

Real 

Degrees 

Line  37 

Littoral  Current 

Real 

Knots 

Line  38 

Number  of  Surf  Lines 

Real 

Line  39 

Surf  Zone  Width 

Real 

Feet 

Line  40 

Wind  Speed 

Real 

Knots 

Line  41 

Wind  Direction 

Real 

Degrees 

Line  42 

Blank  Line 

Character 

Line  43 

Modified  Surf  Index 

Real 

Line  44 

Blank  Line 

Character 

Line  45 

Text  Heading  -  Detailed  Surf  Output 

Character 

Line  46 

Blank  Line 

Character 

Line  47 

Text  Heading  Line 

Character 

Line  48 

Text  Heading  Line 

Character 

Line  49 

Text  Heading  Line  -  Units 

Character 

Line  50 

Blank  Line 

Character 

Line  51 

-EOF  Index  Number 

Integer 

Distance  Offshore 

Real 

Feet 

Water  Depth 

Real 

Feet 

Significant  Breaker  Height 

Real 

Feet 

Maximum  Breaker  Height 

Real 

Feet 

Percent  Breaking  Waves 

Real 

Percent 

Wave  Length 

Real 

Feet 

Littoral  Current 

Real 

Knots 

Page  195 


6.2.2  SURF  3.1  Summary  Output  File 


The  SURF  Summary  Output  File  is  in  the  same  format  as  the  Detailed  Surf  Output  file  in  the 
preceeding  section  without  the  Detailed  output  at  the  end  of  the  file.  The  filename  generated  is 
“*.out”,  where  the  is  replaced  with  the  prefix  of  the  input  file  name. 


Line 

Description 

Type 

Units 

Line  1 

Surf  Header 

Character 

Line  2 

Blank  Line 

Character 

Line  3 

SURF  Model  Version 

Character 

Line  4 

Date  and  Time  of  Forecast 

Character 

Line  5 

Output  File  Name  Information 

Character 

Line  6 

Straight  Coast  Wave  Refraction  Option 

Character 

Line  7 

Landing  Zone  Name 

Character 

Line  8 

Sight  Line 

Real 

Degrees 

Line  9 

Interval 

Real 

Feet 

Line  10 

Starting  Depth 

Real 

Feet 

Line  1 1 

Depth  Profile  Name  or  Beach  Slope 

Character 

Line  12 

Spectrum  Usage  Text 

Or 

Character 

Sea  Wave  Height 

Real 

Feet 

Sea  Period 

Real 

Seconds 

Sea  Direction 

Real 

Degrees 

Line  13 

Spectrum  File  Name 

Or 

Character 

Swell  Wave  Height 

Real 

Feet 

Swell  Period 

Real 

Seconds 

Swell  Direction 

Real 

Degrees 

Line  14 

Wind  Speed 

Real 

Knots 

Line  15 

Wind  Direction 

Real 

Degrees 

Line  16 

Tide  Level 

Real 

Feet 

Line  17 

Blank  Line 

Character 

Line  18 

Internal  Grid  Spacing 

Real 

Feet 

Line  19 

Significant  Wave  Height  from  Input  File 

Real 

Feet 

Line  20 

Significant  Wave  Height  from  Straight  Coast 

Real 

Feet 

Line  21 

Input  Spectrum  Type 

Integer 

— 

Line  22 

Significant  Wave  Height  Offshore 

Real 

Feet 

Line  23 

Stress  Significant  Wave  Height 

Real 

Feet 

Line  24 

Stress  Peak  Frequency 

Real 

Hertz 

Line  25 

Stress  Zero-Crossing  Frequency 

Real 

Hertz 

Line  26 

Stress  Peak  Period 

Real 

Seconds 
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Line 

Description 

Type 

Units 

Line  27 

Percentage  Breaking  Waves  at  Starting  Depth 

Real 

Percent 

Line  28 

Self  Starting  Option 

Character 

Line  29 

Blank  Line 

Character 

Line  30 

Text  Heading  -  Surf  Forecast 

Character 

Line  3 1 

Significant  Breaker  Height 

Real 

Feet 

Line  32 

Maximum  Breaker  Height 

Real 

Feet 

Line  33 

Dominant  Breaker  Period 

Real 

Seconds 

Line  34 

Dominant  Breaker  Type 

Character 

Line  35 

Breaker  Percentages 

Character 

Percent 

Line  36 

Breaker  Angle 

Real 

Degrees 

Line  37 

Littoral  Current 

Real 

Knots 

Line  38 

Number  of  Surf  Lines 

Real 

Line  39 

Surf  Zone  Width 

Real 

Feet 

Line  40 

Wind  Speed 

Real 

Knots 

Line  41 

Wind  Direction 

Real 

Degrees 

Line  42 

Blank  Line 

Character 

Line  43 

Modified  Surf  Index 

Real 
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6.2.3  SURF  3.1  Data  Only  Output  File 

The  data  only  output  file  contains  columnar  data  most  often  used  for  plotting  purposes.  This 
file  was  created  to  ease  the  I/O  reading  for  visual  representation  of  the  values.  This  matrix  of  values 
represents  the  cross-shore  distributions  of  the  variables  defined  in  each  column.  The  filename 
generated  is  “*.dat”,  where  the  is  replaced  with  the  prefix  of  the  input  file  name. 


Description 

Type 

Units 

EOF  Index  Number 

Integer 

Distance  Offshore 

Real 

Feet 

Water  Depth 

Real 

Feet 

Significant  Breaker  Height 

Real 

Feet 

Maximum  Breaker  Height 

Real 

Feet 

Percent  Breaking  Waves 

Real 

Percent 

Wave  Length 

Real 

Feet 

Littoral  Current 

Real 

Knots 
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6.2.4  SURF  3.1  Shallow  Water  Directional  Wave  Spectrum 

The  file  is  only  created  when  Line  10  of  the  Surf  input  file  contains  a  negative  number.  The 
format  of  this  ASCII  text  file  is  a  simple  matrix  of  rows  and  columns.  This  file  has  the  same  file  name 
as  the  output  file  except  that  the  file  extension  will  be  *.dws.  The  first  row  contains  the  center 
frequency  bin  definitions  and  the  first  column  defines  the  wave  direction  bins.  The  heart  of  the  matrix 
is  the  spectral  wave  energy  per  frequency  and  direction  with  the  units  m2/Hz*radians.  This  spectrum 
has  the  same  units  as  the  input  directional  wave  spectrum. 


Line 

Description 

Type 

Units 

Range 

Line  1 

Frequency  Bins 

Real 

Hertz 

0-0.5 

Line  2-EOF 

Wave  Direction,  Wave  Energy 

Real 

Degrees, 

m2/Hz*rad 

0-359 

0-999 
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6.3  Error  Message  Description 


Error  Message 

Subroutine 

Generating 

Error 

Error  105  -  All  input  depths  are  less  than 
starting  depth.  Check  inputs.  Program 
stopped. 

CINDEP 

Error  1 15  -  Opening  Directional  Wave 
Spectrum  File. 

READSPEC 

Error  120  -  Opening  input  file. 

SRFSETUP 

Error  125  -  Opening  of  Input  Depth  File. 

C_IN_DEP 

Error  130  -  Opening  Refraction  File. 

READRFRC 

Error  145  -  Input  depth  profile  has  more  data 
points  than  allowed.  Check  depth  profile. 
Program  stopped. 

CINDEP 

Error  150  -  Large  Internal  grid  spacing. 

Check  depth  profile. 

DEPDRVR 

Error  160  -  No  Convergence. 

NONLIN2 

Error  165  -  No  sediment  type  selected  for 
Equilibrium  Profile. 

EQUILPRF 

Error  170  -  No  Surf. 

SURF 

Suggested  Solution  to  Resolve 
Error 


Decrease  the  starting  depth  in  the 
input  file  -line  9  or  extend  the 
depth  input  profile  farther  offshore. 


Check  Wave  Spectrum  name  in  the 
input  file-  line  5.  Verify  the 
location  of  the  spectrum  file  is  the 
same  as  the  input  file. 


Check  the  name  of  the  input  file 
typed  at  the  command  prompt 
(Surf3.1  <  InputFile)  or  the  name 
typed  during  execution  (Enter  Input 
File  Name). 


Check  Depth  Profile  name  in  the 
input  file  -  line  4.  Verify  the 
location  of  the  depth  file  is  the 
same  as  the  input  file. 


Check  Refraction  name  in  the  input 
file -line  6.  Verify  the  location  of 
the  refraction  file  is  the  same  as  the 
input  file. 


The  maximum  number  of  depth 
points  allowed  is  500.  Modify 
depth  input  file  to  contain  only  500 
depth  values. 


Use  the  Self  Start  Option  in  the 
input  file  -  line  9.  Refer  to  the  Self 
Start  Option  in  Section  6.1.5. 


Smooth  the  input  depth  profile. 


A  Slope/Sediment  Type  was  not  set 
correctly  in  the  input  file  line  8. 

The  value  must  be  inclusive  of 

1  - 10. 
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Check  the  heading  toward  the 
beach  in  the  input  file,  line  7  and 
the  Spectrum  Input  File.  Also, 
there  may  just  be  no  surf  in  the 
area. 


Error  180  -  Problem  gridding  to  output  file. 
Program  stops. 


PRT  OUT 1  Check  that  the  input  depth  profile 

PRT  OUT2  extends  to  the  beach  shoreline  and 

that  the  tide  level  -  line  12  is  not 
too  high. 


Error  185  -  Problem  with  wave  height  values. 


NEW  BRK 


Check  the  input  depth  profile.  The 
data  may  need  to  be  smoothed  due 
to  unusual  slopes.  (Hint:  too  many 
negative  slopes.) 


Error  195  -  Significant  wave  height  outside 
surf  zone  less  than  0.5  ft  -  no  further 
calculations. 

Error  200  -  Surf  forecasts  are  for  situations 
when  waves  are  more  important  than  winds. 
This  is  not  the  case  for  input  waves  and 
winds.  Forecasts  may  not  be  valid. 


S  NOSURF 


Check  the  heading  toward  the 
beach  in  the  input  file  -  line  7. 


S  COEFF 


Check  the  input  wave  and  wind 
conditions  in  the  input  file  -  line  1 1 
and  line  12. 


Error  205  -  Water  edge  not  found.  Check  tide 
and/or  depths.  Program  stopped. 


S  TIDE 


The  input  depth  profile  must  extend 
to  the  beach  including  the  addition 
of  a  tide,  if  specified.  There  must 
be  a  depth  at  either  0.0,  an  onshore 
value,  or  an  elevation. 


Error  210  -  Wave  direction  not  toward  the 
beach  -  no  further  calculations. 


RAD_ST2 


Check  the  heading  toward  the 
beach  in  the  input  file,  line  7  and/or 
the  directional  wave  spectrum  file. 


Error  215  -  Wave  induced  set-up  not 
converging  to  tolerance. 


SETUP 


The  input  depth  profile  must  be 
smoothed. 


Error  220  -  Wave  induced  Set-up  is  not 
converging.  Ending  program. 


MAIN  WAV 


The  input  depth  profile  must  be 
smoothed. 
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6.4  Flowchart  Symbol  Index 


Terminus  (Start,  Return,  End) 


Process  or  Calculation 


Call  to  Subroutine  or  Function 


Connector 


Keyboard  Input 


Loop  Limitation 


Preparation  Prior  to  a  Process  or 
Initialization  of  Variables 


I  Symbol  Connector 


yes  Yes  Connector 

I 


— No- 


No  Connector 
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6.5  Acronyms 


CNMOC 

Commander,  Naval  Meteorology  and  Oceanography  Command 

CSCI 

Computer  Software  Configuration  Item 

csu 

Computer  Software  Unit 

DWS 

Directional  Wave  Spectrum 

EOF 

End  of  File 

Hz 

Hertz 

LHS 

Left  Hand  Side  of  Energy  Balance  Equation 

m 

Meter 

N 

Newton 

MSI 

Modified  Surf  Index 

NRL 

Naval  Research  Laboratory 

OAML 

Oceanographic  and  Atmospheric  Master  Library 

ONR 

Office  of  Naval  Research 

RHS 

Right  Hand  Side  of  Energy  Balance  Equation 

RSM 

Refraction/Shoaling  Matrix 

SPAWAR 

Space  and  Naval  Warfare  Systems  Command 
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